Cosmological Perturbation Theory With 
Background Anisotropic Curvature 

Tom Zlosnik 

E-mail: tzlosnikOperimeterinstitute . ca 

Perimeter Institute for Theoretical Physics, 31 Caroline Street North, Waterloo, 
O '. Ontario N2L 2Y5, Canada 

(N 

3 



Abstract. The theory of cosmological perturbations is extended to spacetimes 
displaying isotropic expansion but anisotropic curvature. The perturbed Einstein 
equation and Boltzmann equations for massless and massive particles are derived 
Ch ' in a general gauge and a decomposition of perturbations into harmonic modes and 

moments is proposed. Generalization to the case where anisotropic expansion is 
also present in the background is discussed. 

cr: 

__I' 1. Introduction 

The observed approximate liomogeneity and isotropy of tlie cosmic microwave 
^ I background (CMB) remains one of the most striking discoveries in cosmology. This 

approximation holds so well that deviations are adequately characterized by a small 
perturbation to the temperature Q = AT/T as well a polarization of the photon field, 



oo 
en 

t^ ■ expected to be of comparable or smaller order on various scales [1,. Both of these 

1^ I perturbational quantities are time dependent, inhomogeneous (having a nontrivial 

f— ^ ■ spatial dependence) and anisotropic (having a nontrivial dependence upon the angle 

that an observer see an incoming photon) . The influx of data from precision cosmology 
has allowed for a more searching question to be asked - that of whether these quantities 
are statistically isotropic. For the case of the temperature perturbation 0, fluctuations 
, , ■ are statistically isotropic if the following identity holds: 

rS 



6(n;i,)e(n;„)) = /(n;„ii,j„) (1) 

i.e. the above expectation value for the product of temperature fluctuations at two 
different directions in the sky n^N and nL^ is a function of only their projection along 
one another, the implication being that there is no other spatial direction present in 
the physics of the generation and evolution of temperature perturbations. The degree 
to which statistical isotropy of CMB temperature fluctuations is true given the data 
has been the focus of considerable debate. There exist indications that observations 
are indeed implying that ([T|) is not exact in our universe on large scales, notably there 
appear to be a number of anomalies that are relatable to the presence of directions 
other than nJ^^ and nLj (see for instance [2] , [2] , [1] ) . 

Given these indications one may speculate as to whether the presence of an 
extra direction in a given CMB anomaly is as a consequence of the presence of fields 
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which impose an actual preferred spatial direction in the universe. This would not 
be unprecedented in cosmology. For instance, the existence of a field to generate a 
primordial period of inflation is conventionally taken to explain the existence of a 
preferred set of clocks on large scales in cosmology [5] . 

If there does existed a preferred spatial direction, this is expected to be encoded 
in the gravitational field either at the background or the perturbed level. Consider 
the following class of background metrics: 

ds2 ^ „2(^) I _^^2 + ( ^ j dz^+ -fcddx^dx'^ j (2) 

led = T^^cd (3) 

The interpretation of the tensor jcd is as follows: 



• 



• 



For the case K > 0, jab is the metric on a 2-sphere with two dimensional Ricci 
scalar ^R — 2 (i.e unit radius) 

For the case K < 0, •jab is the metric on the upper sheet of a 2-hyperboloid of 
two sheets with two dimensional Ricci scalar ^R = —2. 

For the case K = 0, the tensor {l/lKD'^ab is a two dimensional, flat Euclidean 
metric (written in cylindrical coordinates). 



For non-zero values of a{t) and b{t), one may intuitively think of these spaces as 
an infinite spatial line (the 'z direction') with a two-dimensional surface at each point 
along this line- these surfaces may either a flat 2-space (the case K — 0), a 2-sphere 
(the case K > 0), or the upper sheet of a two dimensional hyperboloid of two sheets 
[k <0). 

When ^ 7^ ^ the expansion rate is anisotropic and hence there exists shear in 
the congruence of timelike geodesies labelled by proper-time t = J adt. The form 
of ([2]) additionally implies that this congruence has vanishing vorticity. Cosmological 
perturbation theory for the gravitational field coupled to a scalar field for the case 
with shear and K — has been developed in detail ([B],[7],[H])- The case where shear 
is present and K y^ has also been considered at the level of the background [9] . The 
authors explicitly consider a scenario of primordial anisotropic inflation, followed by 
isotropic, flat FRW expansion- such expansion corresponds to the choice K = 0,a = b 
in ([2]). The anisotropy is then encoded in the perturbations and not the background. 
However, we will see in Section IH] that this scenario is also straightforwardly addressed 
in choosing a '2+1' decomposition of the background space as a result of how this 
allowance for spatial anisotropy is encoded into the perturbation theory. 

Alternatively, deviations from Friedmann-Robertson- Walker (FRW) geometry 
may occur in the absence of shear but in the presence of anisotropic curvature (i.e. 
\K\ y^ 0). This case was discussed |I0] and further developed in [ll],[12]. Crucially, 
in the absence of shear the metric ^ is conformally static and therefore admits 
a homogeneous and isotropic solution for the CMB at the background level |13) . 
Anisotropic cosmological models which lack the conformally static property typically 
have fine tuning problems as they they induce anisotropics in the CMB radiation which 
are known to be consistently regarded as small perturbations around a homogeneous, 
isotropic solution. Clearly the models are also immune to constraints on directionality 
of the expansion rate so in the shearless case departures for the predictions of the 
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FRW model begin instead via the anisotropy of the luminosity distance relation which 
can used to constrain K by comparing with data from supernovae. It was found [T3] 
that constraints on the quantity fix (i-e the contribution to the critical density of 
components with energy density scaling as a~^) from the Union SnE dataset |15) are 
comparable to those obtained in the FRW model. 

The intention of this paper is to develop the theory of cosmological perturbations 
in the case where the cosmological background belongs to the class ([2]), is shear 
free but may posses anisotropic curvature. In Section [2] we discuss the harmonic 
decomposition of perturbations and present the perturbed Einstein tensor in a general 
gauge. Following this we consider perturbations in the matter sectors. Firstly in 
Section [3] we consider perturbations to a putative field supporting the background 
shear-free condition. Then in Section 2] we consider the Boltzmann equation for 
massless particles, and in doing so propose a decomposition of the angular dependence 
of the particle distribution function in terms of moments. In Section[S]we compare the 
proposed moment expansion to that conventionally employed in the FRW limit, and in 
Section [6] we discuss the consequences that the mere presence of anisotropic curvature 
at the background level may have on the evolution of of the photon perturbation. 
Then, in Section [71 we consider the Boltzmann equation in the case of massive 
particles. In Section [8] we present expressions for the angular correlation function for 
polar temperature perturbations and in Section [3] uses these results for comparison to 
previous results in the literature. The extension of the formalism to the case where 
shear is present in the background is discussed in Section [10] and conclusions are 
presented in Section [TTl 

Additionally there are a number of appendices: some necessary mathematical 
background to the calculations is presented in [Appendix A| the explicit calculation of 



the polar angular correlation function is given in Appendix B| the form of spacetime 
gauge transformations and the construction of gauge invariant perturbations is given 
in [Appendix C[ comparison to scalar- vector-tensor perturbations in the FRW limit is 



given in Appendix D[ and a table of notation is given in [Appendix E[ 



2. Cosmological perturbations 

For the shearless case, as in the case of linear perturbations around an FRW spacetime, 
the problem of deducing the spacetime metric, if sufficiently 'close' to the background 

metric, is equivalent to a problem of examining the evolution of a set of fields on a 

(c) 
'co-moving spacetime', this spacetime defined by the metric 5/1^ : 

gl^Jdx^dx'' = -dt^+ dz^ + -fcddx'^dx'^ (4) 

In the case of linear perturbations around an FRW background, surfaces of 
constant conformal time on the comoving spacetime are homogeneous and isotropic. 
This isotropy allows for a simple decomposition of fields into scalar, vector, and tensor 
perturbations, which evolve independently of one another |16| . This decomposition 
is not appropriate in the case ([4]) where K ^ 0, due to the inherent anisotropy of 
surfaces of constant conformal time. If attempting to do so, one would find coupling 
between scalar, vector, and tensor modes, so considerably complicating the analysis 

Our approach is to follow ([IH],[I1]) and take advantage of a more limited isotropy 
present in the spacetime - i.e. the isotropy of surfaces of constant (i, z). Due to this 



Cosmological Perturbation Theory With Background Anisotropic Curvature 4 

isotropy, we may decompose a function j/(i, z, x") existing on the co-moving spacetime 
as follows: 

j/(i,z,a:«)=^y(i,z)f^(a:«) (5) 

k.rn 

Unless otherwise stated, coordinate indices a, 6, c. refer exclusively to components 
in the 2-surface. The functions £^{x°') form a complete set of eigenfunctions of the 
Laplace-Beltrami operator on surfaces of constant (i, z), and indeed for a general value 
of the curvature K, we define the complex label k via the following equation: 

V,V''f^,(x'') = ~kk*£!;,ix'^) (6) 

where Vq is the derivative operator compatible with the metric jcd- 
Due to our choice of coordinates, the label m is always discrete, though with 
a range dependent upon the curvature of the surface, whilst k may be continuous 
or discrete depending upon the curvature of the surface In the continuous case the 
summation over k must be replaced by an integration with appropriate measure. For 
instance, for K > Q, the functions are the spherical harmonics l^Jj, and the eigenvalues 
kk* are equal to |iir|^(/ + 1). See Appendix... for discussion of the corresponding 
functions in the open and flat case. 

Similarly, a vector field ya{t, z,x^) (i.e. a two component vector field existing on 
surfaces of constant (t, z)) may be decomposed as follows: 

ya{t,z,x') = Yl (y^'''^Ht,2)Vaf^(x'') + y(^^)(t,z)V„f^(a;^)) (7) 

k,m 

Where we have now introduced the differential operator operator Vq = eJ'Wb 
where eab is the volume form on the co-moving 2-surface. The labels P and A refer to 
polar and axial terms respectively. It may be checked that 'VaSmi^'') ^-nd Wa£m{^^) 
are eigenvectors of the Laplace-Beltrami operator on the comoving 2-surface. 

We shall find that, much as scalar, vector, and tensor modes evolve independently 
of one another in the FRW case, the polar and axial modes evolve independently for 
small perturbations around the metric ([2]). 

Finally a symmetric tensor field yabit, z^x"^) existing on the 2-surface may be 
decomposed as follows: 

yab{t,z,x^) = ^(y(^^i)(z,t)7afcf^ + y(^^'H;^,t)VaVfcf^ 

+ y(^^)(z,t)V(,Vfc)£^) (8) 

Where P and A again denote polar and axial modes, and V a^b^m ^'^d VaV^f^ 
are eigentensors of the Laplace-Beltrami operator on the co-moving 2-surface. Note 
that the appropriate range of k need not be the same in the respective decomposition 
of functions, vectors, and tensors on the 2-surface. 

The decompositions ([5]),(l7]), and ^ then allow for the decomposition of a general 
perturbation 5g^i, to the metric ([2|) into polar and axial components as follows: 
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%. = Sgi!:^ + Sgii^ (9) 



where 



xz^ 



Sgl^Ux^'dx" = a^ ^{ViyZ, t)edt^ + E{z, t)£dzdt + F{z, t)8dz 

A;, 771 

+ B{z, t)Vd£dzdx'^ + C{z, t)Wd£dtdx'^ 

+ U{z, t)-iab£dx°'dx^ + X{z, t)VaVb£dx°'dx^) (10) 

and 

Sgl'^Ux^dx" = a^ ^{R{z , t)W a£ dtdx" + S{z,t)V aSdzdx" 

+ 2Q{z,t)V(^a^i,)£dx°'dx^) (11) 

where the labels k and m on 5,^^ have been suppressed for notational compactness. 

Using the expansions p^ and ([TT|) along with the eigenvalue equation, we will 
find that the Einstein equations may be expanded in terms of sums of eigenfunctions, 
eigenvectors, and eigentensors, with coefficients of which being differential equations 
for the harmonic modes of the expansion. These differential equations will contain 
partial derivatives with respect to time and the z coordinate. It is desirable to 
transform this system into a set of ordinary derivatives with respect to time by instead 
looking at the evolution of the Fourier transform of the function yz,t^ defined as follows: 

y{z,t) = ^j y^{t)e'^^dw (12) 

2.1. Polar perturbed Einstein tensor 

We now detail the components of the polar components of the perturbed Einstein 
tensor 5G in terms of these Fourier modes. We shall use a prime to denote 

ordinary derivatives with respect to conformal time t and in our expressions the the 
appropriate multiplication by e™^ and sum/integral over labels k,m,w is implicit. 
For notational compactness we denote the real number kk* by k^. 

6G'-^^\ = -^(2niwE + Sn^V - 2nk^C + iwk^B 
^ a-' 

- -k^F - -HE' - w^U - 2nU' - -k^U + KU 

+ -k'^w^X + nk'^X')£ (13) 



a^' 2 



5G^^^\ = -i^{-^{E - iwC -B' + iwX') - UiwV + iwU')£ (14) 
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5G^'"'\ = ^iJiV + 1{-2H' + 4- - k')V£ 
^ a^ 2 a 



- k^{C' + 2nC)e - ([/" + 2nU' -KU+ -k^U)£ 

+ k^{^X" + nX'))£ (15) 



^G^^)* ^J-(-iwE^2'HV+(-w^C + 2KC) 
" 2a^ 

- iwB' + F' + U' - KX')Va£ (16) 



+ {B" + 2HB' - 2KB) - iwU + iwKX)Va£ (17) 



^G^'^^\ = A ((^" + 2H^')7V + F{V-V,£ + e-f\£)) 



2a2 



— r((-2-HV + 4—1^ - w^V + 2-HT/')7V 
2a^ a 



V(k''-f\£ + VVbS)) + -^iwB (V^Vbf + k^-/\£) 

4j(C' + 2nC) (V^Vfcf + fc27"tf ) + \iw {2HE + E') -i\E 

±^{^w^U-U"-2nU')j\£ 

^ {X" + w^X + 2HX) (V' Vfcf + k'-i\S) (18) 



The perturbation (llSp has two tensor components: a sum of perturbations 
proportional to ^""^E and a sum proportional to V°Vbf . One may isolate the latter 

contraction with traceless components. 



components by operating on (J18l) with VqV'' — i7'^[,V2 which is nonzero only after 



2.2. Axial Perturbed Einstein Tensor 

We now consider the axial perturbations to the Einstein tensor: 

5G^''\ = ^(^^^^^i? + K?, - ZS' + ?Q' - \kQ')V^E (19) 



5G^^^\ - ^(-^(^' + 2^^) + \s" + ns' + ^S- l-KS 

a"^ 2 2 2 2 

-^Q+yO^)V.f (20) 

Observe that there are axial and polar contributions to SG^^. However, their 
decoupling is guaranteed by orthogonality of vectors \/aS and Va£- 
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1,1,1 , iw 

--{R' + 2nR) + — 



6G^^^\ = ^{-^{R' + 2HR) + ^S 



+ iQ" + u?Q + TH^Q - A—Q + 2?^Q'))7°'V(,Vm (21) 

a 

Again, orthogonality of axial and polar contributions to (^G"^ follows from the 
7"'^V(cVfc) having vanishing contraction with ^J'£ and VaV''.?. 

Given the perturbed Einstein tensor it is necessary then to obtain the perturbed 
stress energy tensors of matter fields in the universe, as well as their owned perturbed 
evolutions. We will obtain both of these things by a treatment of the Boltzmann 
equation for massless and massive particles in Sections U and [T] In the following 
section we discuss the perturbations of an explicit example of a field whose presence 
can allow for a shearless example of the spacetime ([2]) to exist as a consistent solution 
to the Einstein and matter field equations. 

3. Scalar field 

As noted in [TU], a shearless example of © with K ^ Q cannot be be a solution of 
the Einstein equations sourced by a perfect fluid. Rather, it is necessary for there 
to exist a matter source in the Einstein equations with an anisotropic stress at the 
background level. Indeed it may be shown that if the matter source may be described 
by the following fluid stress energy tensor: 

V = pU^U, + Ph^, + LV^V, (22) 

where U^, is the fluid 4- velocity, h^^, is the metric on surfaces of constant cosmic 
time, and V^ is a unit vector (with respect to ([5])) along the {8^)^ direction. In turn, 
the anisotropic pressure L must satisfy the following relation in order for ([5]) to be a 
solution: 

K 

L = -- (23) 

One would not expect there to be a unique field/particle origin of a matter 
candidate that could satisfy (1^ . but it is instructive to consider an explicit realization 
of the above condition. This was provided in |14) via a massless canonical two-form 
field. Equivalently, such a theory is equivalent to a massless scalar field with the 
following action [20] : 



C 



S = -— — / d^x^/q^V'>V^0 (24) 



167rG 



We consider the following ansatz for the scalar field 

(l){z,x'',t)^ — + W{z,x'',t) (25) 

where the field M^ is a small perturbation. The ansatz for the background field 
configuration guarantees that the background metric takes the shearless, curved form 
if zq is related to C and K as follows: 
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-2--K (26) 

Up to first order in perturbations, the scalar field stress energy tensor T'^^, has 
the following nonvanishing components, up to first order in perturbations: 

where the summation is understood here to denote a sum/integral over the 
appropriate harmonic labels k,m,w 

T* = -^Y ( — - w] Se'""' (29) 



J2 WVaSe'""' (30) 



" a^zo 



Additionally, the scalar field equation of motion is given by 

,^,,, ^, . , 9 ,9,^,, iwU iwF iwk^X E' 

W" + 2nW' + {u? + k^)W = + — 32) 

zo 2zo 2zo zq 

Note that the scalar field is coupled solely to polar perturbations in its own field 
equation and in its contribution to the Einstein equations. The background energy 
density scales as a curvature component would. The restrictions on allowable values 
of K (and thus C via ([26l) due to probes of the background spacetime imply that 
as a component of the cosmological energy density the scalar field will only become 
non- negligible at late times [H]. Note that the perturbed stress energy tensor also 
scales the density of a curvature component multiplied by a small perturbation. 

4. Boltzmann Equation for massless particles 

We expect that in a realistic universe it is appropriate to describe a field ^ defined 
on the spacetime @ in terms of its particle content via a quantity /*(x^, P^^) called 
the distribution function, which refiects particle number density in phase space at a 
time t. This function is expected to obey the Boltzmann equation, which is used to 
calculate the function's total time derivative: 

dt dt dx' dt dP^ dt ^ ' 
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where i denotes the collective labels z and coordinates on the co-moving 2-surface, 
and the partial derivative is evaluated at constant phase-space point. The term C is a 
coUisional term. 

We initially consider the collisionless (C — 0) limit of the Boltzmann equation for 
a massless field, described by a distribution function / and restrict ourselves to the 
case where the massless field is the photon field (though the results of course will apply 
to any massless field). A considerable simplification that follows from the condition 
C = is that the particles will follow geodesies of the perturbed spacetime. The 
dependence of a function upon z and x"' is familiar from the preceding sections, but 
the dependence on 4-momcntum P^ is new. For a massless particle, this momentum 
satisfies a null constraint: 



g^^P^^P" = (34) 

Therefore, only three components of the 4- vector P^ may vary independently. We 
shall find it convenient throughout to use the constraint (I34[) to eliminate the time- 
component P* of the photon 4-inomentum. We write the 4-momentum as follows: 

P^d^^-{P'dt+pP'd,) (35) 

a 

P* = aP* (36) 

P'di =ad,+pVl-a^p''da (37) 

What are the interpretations oi p,a and p"? Calculating the norm of the spatial 
according to the background spatial metric hi 
moving background spatial 3-metric, we have that: 



part according to the background spatial metric hij = a'^hl, , where h], is the co- 



h,,P'P^ = h^^p^P'P^ - p^a^ + p'(l - a')7a6pV (38) 

(C) 

where /i^ is the co-moving background spatial metric. If we identify p"' as 
components of a unit vector living in the background co-moving 2-surface i.e satisfying 
labP°'P^ = 1 then clearly p^ = hijP^P^ i.e. it the norm-squared of the physical spatial 
momentum of the photon according to the background spatial metric. The variable 
a is then interpreted as the projection of the co-moving spatial momentum alongside 
a vector in the z direction which is unit with respect to the co-moving background 
3-metric. The value of using the variables a and p° is that they more clearly reflect 
the underlying symmetry of the background spacetime. 

The total time derivative df /dt is then decomposed as follows: 

dt dt dx' dt dp dt da dt dip") dt ^ ' 

It is useful first to express the derivative ^ in terms of the photon momentum. 
Firstly we note that we may define the particle four-momentum in terms of the affine 
parameter of null trajectories A as follows P^ = dx^ jdX Therefore we have that: 

dx^ jdt = {dx^ /d\)/{dx^/dX) ^ P^ /P^ = pP^ /P* (40) 
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Next we eliminate P* via the null-constraint p4p . Up to first order in metric 
perturbations this equation reads: 



- - (1 + V){P*f + 2p{aE + {Ba + Ra)fW^~ol^)P' 

+ p^ + f/(l - a^)p^ + FaV + 2(5a + ha)p''a^/l ~ a^p"^ 

+ {l-a')p\Xab + 2Qab)rp'' (41) 

Where for notational compactness the metric perturbations here are the actual 
metric perturbations i.e. the appropriate sum, integral of all the harmonic 
components/Fourier modes of the perturbation. Solving this equation for P* we have 
up to linear order in perturbations that: 

P* - p+p{aE + [Ca + Pa)pVl-a2 + ^(l - O?) 

F 



-a^ + iSa + BaWa^l-a-' 

+ i(l-a2)(x,, + 2Qa6)pV-y) (42) 

Consequently, using (|4]) in conjunction with ()40p . we may express dx^ /dt entirely 
in terms of the variables p,a, p*", and metric perturbations. 

Another term that we must calculate is dp/dt. From the time component of the 
geodesic equation we have that: 

dP* P* d /P*\ , „ „ 

We may use this equation in conjunction with the solution and the observation 
that the geodesic equations (as may be checked) imply that da/dt and dp"^ /dt are zero 
at the background level to find the following expression for dp/dt: 

j-= -pH-J2 ^ia^9,F + 2aE' + ad,V + 2(1 - a^)U' 
+ 2a^F' + a{l-a'^)d^U)£ 
VT 



J2^^^^iia^-l)U-2a'd,B 
- 2C" - AaB' -V - a^F)p''Va£ 
+ ^(1 - a^)p{-Ba -X' - ]-d,Xa)rp''yaVb£ 



J2 o "' (1 - «')p^rpy V,VbV,£ 



+ J2 Vl - a^p{-R' - 2aS' - a^d,S)p''\Ia£ 

+ ^(1 - a')p{-ad,Q - 2Q' - aS)p''p''V aVbS 

- Y, Vl-a2(i _ a^)pQp-pbf^^Vf,Vc£ (44) 

where the sum denotes the appropriate sums/integrals over harmonic components 
and the plane wave e™^ is implicit. 
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We now decompose / into an unperturbed component fo{p,t) and a perturbed 
component Sf{Pi,x^,t). Note that we have assumed that the perturbation has a 
neghgible dependence upon the photon momentum magnitude. 

In a reahstic universe, the unperturbed distribution function /o is expected to be 
of the Bose-Einstein form. Motivated by this, we introduce a new field Q(Pi,x'^,t) of 
first order in smallness as follows: 

/ = P^ (45) 

gT(t)(i+e) _ I 

where T is the background field temperature and we recall that p is the magnitude 
of the physical photon S-monicntum. From the background geodesic equation we 
recover the solution p = pa/a where po is a constant. A simple calculation shows that 
the background Boltzmann equation df/dt = is satisfied if T = To/a, a familiar 
result. 

Our definition of the fields 8 and Sf implies that: 

Sf = -P^Q (46) 

We now look to find the collisionless Boltzmann equation to first order in 
perturbations. We note that as each of ^, ^, gfL) , and -^-^ are quantities of 
order 1 or greater in perturbations, the combined terms appearing in (|39p may be 
neglected to first order. We have then that: 

dt \ dp J dt dp 

where '°' and —rr- denote background and perturbed time derivatives of p, 
respectively. Furthermore, it may be checked that 

dfo d ( dfo\ rfp(o) ..„. 

-""^'^-d-pV^j^r ('') 

The ensuing cancellation of terms implies that: 

Additionally, dx^ /dt = pP^/P* = P* to zeroth order. Then, reading off -^^ from 
(l44l) we have: 



^-{a^iwF + 2c 
2a^F' + ail~a^)iwU)£ 



Y^ -{a^iwF + 2aE' + aiwV + 2(1 - a^)U' 
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- Y, ^^~"' ((q' - l)U - 2a^iwB 

- 2C' - AaB' -V - a^F)P''Va£ 

~ ^(1 - a^){-Ba -X' ~ iz«;Xa)py VaVftf 



- Y^ \/l-a^{-R' - 2aS' - a^iwS)p''V a£ 

- ^(1 - a^){-aiwQ - 2Q' - aS)p"' p^V aV b£ 

+ ^ v/l-a2(l - a2)QpypCVaVbVc£' (50) 

Much as it is often more convenient to decompose a position dependent field 
into harmonic components, it may be additionally convenient to decompose a field's 
dependence upon the unit vector P* into 'moments' so that the partial differential 
equation (|50p may be cast as a system of ordinary differential equations in conformal 
time for a set of moments for each harmonic mode. For the field 0, the appropriate 
decomposition of this dependence shall be expected to depend on the symmetries of 
the problem and the nature of the source terms in the equation. 

We first note that the source terms from the metric perturbations are a 
combination of sums of polar and axial harmonic modes. It shall be useful then 
to split Q{x\P\t) into contributions Q^P\x\P\t) and Q'^'^\x\P\t) such that 
@{x\P\t) = Q'^P\x\P\t) + Q'^^\x\P\t). As the labels suggest, these are polar 
and axial contributions to the temperature perturbation. We will now consider these 
contributions in detail. 

4-.1. The polar temperature perturbation 

We propose that for a given curvature K, the field Q^^\x^ , F*, t) may be decomposed 
as follows: 

Q^P\x\P\t)^ Y E 4^AKfc)0MAKfc.i)^MA (51) 

A,0<M<Ak,m,w 



where 



Bi^{{z,w,k,m,a,pn = ' "^^^^"^"^ ^J^^^ ^^"■■"^^"' 



(52) 

Va6..Af = VaVfc....VM (53) 

P^j{a) is an associated Legendre function, A is a non- negative integer and we 
introduce the Mth rank tensor C°'b..M ^ which wc shall call a polar Chebyshev tensor 

(p) 
for reasons made clear below. This choice of basis function B\jj^ is strongly motivated 

by the treatment of cosmological perturbations on curved FRW backgrounds, where 
one decomposes the angular part of the photon distribution (i.e. a and p°) in terms 
of Legendre tensors (see [H] , [H] , [23] ) ■ 

Note that the value k shall refer the positive square root of kk* . The properties 
of this tensor will be discussed below. We shall seek to ultimately relate the moments 
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Q}^'^{w,k,t) to observables, and so the functions c}^'^{w,k) shall be chosen so as to 
simplify the resulting expressions. 

4-1.1. Polar Chebyshev Tensors We define polar Chebyshev tensors via the following 
scalar, vector and recursion relation to generate tensors of higher rank: 

C = 1 (54) 

C = p" (55) 

C-b-N+i ^ 2p('^C''-^) - -^(''^c-^-i) (56) 

where brackets denote symmetrization and the notation C""^'^ denotes a tensor 
with N indices. The tensors satisfy the following properties 

(j.a.b.N ^ (j.b.a.N (57) 

C-^-^'-^'lab = (58) 

C-^-% = C-^-i (59) 

C"^^^-^eaefce,..eAr = CAr(e"pa) = cos(iVcos-i(e"pa)) (60) 

where e" is a unit vector with respect to the comoving ^ab and C'N{e°"Pa) is a 
Chebyshev polynomial of the first kind. 

For notational convenience we define the following function: 

viP = {-k)-'Hl - a^)^C^''-'''Vab..M£'^ (61) 

If Q^^\x'\P'',t) allows the decomposition (|5ip . a relevant quantity to evaluate 
will be the directional derivative y/1 — oP-ff^V a acting on Bj^l^ in ((50)) . The part of 
this derivative independent of a may be developed as follows: 

+ ^ E^'^C'^-^ •'^^a./'V..,..^-if,^„ (62) 

^ar''bc..NY7 ck , ^ar^bcN^l ck 

" ]\l ,l P ^ "^abc.N+lC^ + ^ P O Vbac.N+lC^ 

I ^ ~ ^ ^ar-ibc..N+lY7 ck 

+ ]^^E^-'-""'V-^-^^ (63) 

i 



where Rg^i,t is the Riemann curvature tensor on the comoving 2-surface and we 
use the curvature convention of |24| . After the appropriate number of exchanges of 
indices, we recover: 

(1 - a2)^^a^6c..7V^^^^^^^^^^fe^ ^ (^ _ ^2^E^piacbc..N)y^^^^^^_^^^k^ 

+ ^i^iV(iV - l)K{-k)''-'VM-i 

(64) 
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The first term on the right hand side is closely related to a term from 
the recursion relation ([55)) . We then seek to express the other term from ([55)1 . 
^(a6(^c...Af-i)y^^^ ^^^^ in terms of our mode functions. This is indeed possible, a 
lengthy calculation yielding the following result: 

(65) 
Therefore, collecting results: 

/ 1 r'iab..N+l . (ab(^c..N-l) 

\1 — a^ 



;:; Vaf,..JV+l'^„ 



Vi 



a 



2 



3 



■N{N ~l)K{^k)-'VN 



-1 






yr^^2 



where 



2 ^VT^7?^+i "^^ " l^JV^lJ 



^^^1-7V(7V-1)^ (67) 



(66) 



This implies that 



VT^^r^aB^ - \ «l + -Ui - »'w'.Z) (-i)-l^iM.™. 



= ^nl (a{A -{N- l))4''2i,^ ~{A + N~ l)SJv''-i.A-i) 

^ 2 2A + 1 V^^+i-^-i '^N+uA+i) 

Furthermore from the recurrence relations for associated Legendre functions we 
have that: 

(2A + l)aB\^^\ = iA-N + l)4^]i+i + {A + iV)Si,^Li (69) 

and therefore 
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vT^^p^^^aBn = \ {<i + -^(1 - «^)v^^j,) (-1)"^^^ 



k 1 



2 (2^+1) 






N+1,A+1 



k 2 (A-iV+l)(A-JV + 2) ^(p) 
+ 2^ (2A+1) ^A^-i.A+i 

_ fc (A + iV)(A + iV-l) (P) 

2«A^ (2^+1) ^A'-i.A-i (7U) 

The directional derivative in the Boltzmann equation also involves contributions 
from variations of Bma along the z direction. This is position dependence is more 
simple due to the plane wave decomposition along this direction. Again utilizing 
recurrence relations for associated Legendre functions, we have that: 

"^^^Ka - * (^XTT) ((^ - ^ + ^)^nZ+i + (^ + ^)<Li) (71) 

Collecting terms, we have that the three dimensional directional derivative of the 
basis function Bj^^J^ is given by: 

B=A+1 M=N+1 

j(-P) — Y^ Y^ ^AiN-l-M,A-l-B)/^T j\r(-P) 



where 



B=A-1M=N-1 



-kK% {A+N){A+N-1} „ kK% {A-N+l)(A~N+2} 

2 (2A+1) ^ 2 (2A+1) 

A// = I iw{A+N) p. iw{A~-N+l} 

' (2A+1) " (2A+1) 

k 1 n k 1 

2 2A+1 " 2 2A+1 

We designate rows and columns so that the 00 index position (upper left in the 
matrix) corresponds to Af = A'^ — 1 and B = A — 1. Advancing down a column 
increases M by 1, advancing to the right across a row advances i? by 1. 

Therefore we have that 

p^v,e(^)(x\F\t)= E E E4>^'lKM) 

w,k,mA,AI,B N=0 
._^(JV-l-M,A-l-i3)(^^)g(P)^ (73) 

This is a crucial result as it enables the partial differentials to be replaced in (f50| 
by sums of other moments of 0'^^(a;*,P^i). Additionally, we note that all of the 
polar metric source terms in ([50]) may be expressed in terms of the SJy^ functions. 
Consequently, the polar Boltzmann equation (j50p may be expressed in terms of a 

(p) 
sum of BpfJ^ functions. It is desirable now to find an operator that, when applied 

to the mode and moment expansion of Q^^'{x^,P^,t) yields a term involving only a 

single moment Qj^J^{w, I, t). This will enable us to express the partial time derivative 

of 8*^^-' (a;*, P*, t) the time derivative of the moment 0]y^(w, /, t) along with terms of 

metric perturbations and other moments of the photon distribution. Collectively these 

ordinary differential equations in time will be equivalent to (|50|) . 
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4.1.2. A polar projection operator Consider a field V{x'^,P^,t) may be decomposed 
into polar moments Vj^Aik, ,w,t) as follows: 

V{P\x\t)^ ^ ^y^^(A:,u;,i)SJ^^] (74) 



and recall that P' — {a, y/l — a^p"). 

We seek to find an operator that can operate on V{P^,x'^,t) to isolate a single 
moment VANik,w,t). 

Let us consider the operator O'-^^M' ,M,B,k' ,m' ,w') defined as: 

m- 

(75) 

where d^S is an area element in the co-moving 2-surface and dJ7 = —dadx is a 
solid angle element in the space of photon momentum i.e. % is the angle between p° 
and an arbitrary unit vector at the origin of the momentum space surface orthogonal 
to {dzY. The meaning of the newly introduced parameter tp will be discussed later. 
Consider the following operation: 



O^^^V{P\x\t)= Yl VNAik,w,t)0^^\N,M,B,k',w',m)Bi 



■'na 

A,N,k,m,m\w 



(76) 



We will evaluate this integral in stages. Examining (jllip and ([71)1 . we see that 
5]) will contain the following integral: 



/' 



^^fjab..M(jcd..N (77) 

We introduce a tensor W"-b...m^^ ^^^ defined as follows: 

(jab..m ^ VF'^fc-™ ^ ^^pc^rf..._p™ (78) 

For instance: 

Therefore ([77]) may be written as follows: 

^^C'^b..Mcef..N ^ ^a...Af^^ ^^^e/..A.^^ ^ ^ dxp"" P^ . .p'' p" (80) 

The integral over p°^p^ .. in (j80p may be evaluated by expressing each p" in terms 
of the introduced angle x an a local orthonormal basis i.e. p" = (cosXiSin^)- The 
integral may then be evaluated using standard methods, yielding: 
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P p ...p p ax 2(M+^)(((M + iV)/2)!)2 ^ > 

where we have defined the tensor ^i<^bc..m) ^ ^(ah^cd^ _^.m)^ j-j- jg |;a,ken that the 
result is only nonvanishing when (M + N) is an even number. 

By the definition of the tensor yy'^'>---™^^ ^^ jt must be traceless with respect to 
contraction with 7*^^ across any two lower indices. Therefore, (j77p is only nonvanishing 
if ry(ab..(M+N)) jg contracted with one index of each component "/'^^ acting on one of 
the two W tensors, and the other acting on the other. This implies that M = N. 
By a counting argument, there're 2'^Q\{Q)\ permutations of -y°-''--^Q that fulfil this 
condition. Therefore: 

i^,riab..Mr<wf..M Tjrab..M ^^ref..M (2M). „,,(cfl..2M) 



ff;i..M2(2Af)-i((^/)r 



^ ^y^ab..M.cd..My^,ef..M (2M)! 2^^(M)!(M)! 

cd..M2(2Af)-l((^,/)!)2 (2M)! 

= ^W^'-''^^'-^'W^f-'L.M (82) 

Given the condition that N — M, we may now simply evaluate the integral from 
standard results: 

daPUa)P^Aa) ^ - ^J^J^^^^^^^-b (83) 

Similarly the following integral holds: 

+00 

S{w - w') (84) 

We have then reduced the expression to: 



S{w - w')2 ^^ ^ \ J^^ ^ '■ ^^J abSmnt^m (85) 



{A + M)\ ^ ^ 27r 
■(2A+1)(A-M)!^^^^"^"'2 

The final integral to evaluate is 



This calculation is simplified by the introduction of some new quantities: 
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4-1.3. The A vectors Consider the vectors (ee)" — {dg)"' and (e^)" — 
{^1 \/lH>){(^(j})°' 1 where the coordinates 9 and cj) are angular coordinates on the unit 
sphere for the case K > Q, cyhndrical coordinates for the case K — Q, and hyperbohc 
angular coordinates in the case K > Q. This set of vectors are orthonormal with respect 
to the metric 7°^ except at the origin. We define the following complex vectors: 

(A+r^^i=-((e,r-^(e^n (87) 

(A-r^^^((e,r + z(e^r) (88) 

Note that in the case \K\ = 0, the relations ((ST]) and dHH]) reduce to (A^'^)" = 
(e^)" =F i{^4>Y where r and (j) are cylindrical coordinates on a flat 2-surface. 
The following identities then hold for arbitrary values of K: 

{eoT =^-^^^{iA+r + iA-r) (89) 

i\K\-^^^ 

{e,r =±^^{{A+r~{A-r) (9o) 

(A+)''(A+), = (A-)'^(A-)=0 (91) 

(A-)''(A+), = (A+)"(A-) = 1 (92) 

lab ^{A+)a{A-)t + {A-)a{A+h (93) 

^eab ^ (A+)aiA-)b - {A-)a(A+)b (94) 

Furthermore, for compactness of notation we define tensors 

j^{±)ab..M ^ (A±)'^(A±)^.(A±)^ (95) 

We may use (A )" as basis vectors with which to construct tensors. For instance, 
a component X+^ of a tensor Xab is given by {A~)°'{A~)°'Xab- As a result of 
the properties ([91]) and (|92l) . a tensor Xab.. expressed in terms of (A )° has only 
components which are traceless with respect to any possible contraction of indices 
with 7°'': the components X±±,,. 

Consider a tensor operator which acts on another tensor to produce a tensor 
proportional to part of the tensor which is traceless for all possible contractions and 
symmetric for all exchanges of indices. VV^^J",^ is an example of such a tensor but 
need not be the only possibility- we shall see that a different tensor may be constructed 
using the metric and the volume form on the two surface. The operation of on a 
tensor Xab..M will yield the following result: 

{0r''-^',,..M{XUb..M = P+ ■ (X++..m)(A(+)),,..m 

+ /3_-(X__..Af)(A(-)),d..M (96) 

where /3+ and /3_ are constants independent of X. We now find /3+ and /?_ for 
the case (0)''''-^^cd..M = ^°'''"^cd..M- % its definition we have that: 

W"^''-^',a..MPaPb.-PM = C,d..M (97) 
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For notational convenience we introduce the notation p°^-*^ = p'^p^..p^^ 
By explicit calculation, we have 



Ca^P+iA+)a+P-iA-)a (98) 

Then repeatedly applying the recurrence relation (I56p . we then have that for 
m> 1: 

Ccd..M = 2^'-^p++„Mi+){A + U..M + 2*^-ip__..M(-)(A-)cd..M (99) 
This implies that /3(p)+ = 2*^"\ /3(p)_ = 2^^^^ where M > 1, and hence: 

W"" ■■ cd..M^ab..M£m = 2 ' ^^ ■ "^ ++ ..uSmi-^ )cd..M 

+ 2*^-1 • V__..Mf™(A(-))ed..M (100) 

The integral (I130p . now accompanied by the integration over tp reduces to: 

'(^'''-^'^^dij /d2522*^-i((-i)"'+'"v++-*^rl,v__..Mf-™ 



+ V++-*X'5V__..Mf^) (101) 

where we have used the fact that E*'^ = {—l)P£'lp for each of our cases. 
The total expression is then proportional to: 

'\^ V (2A+l)(A-M)!(-fc)- di: dS 

. ((-i)™'+™v++-'^-^£:il,v^_..Mf^„ + v++-^'£*Jv^^..M£ik) 

.i,™' V (2A+l)(A-M)!(-fc)-^ J ^J 
■ ((-l)-(™'+™) + l)(V++-*^f;„'5'v__..Mf^) (102) 

where we have used the fact that for a value of \p\, the total expression is a sum 
over —\p\ and +|p|, enabling the above rearrangement terms in the summation above. 

We first consider the closed K > case. By the results of [Appendix A[ we have 
that 

,++..A/^5^_ ofe ^ (\K}_Y (^-^^r^) T^'}*T^ - -^-"'+™ 



■ MO„i — \ -;^ I \ — ] J-m',A/-'-m,A/' 



l[ViL + q){L-q + l)\ il[V{L' + q){L'-q+l) 
\q=0 ) \q=a J 
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The matrix elements T_ ,'\^,T^ j^j are orthogonal with respect to integration 
over the group [25] i.e. 

'^-m*M'^-m,M Sm OdlpdfjydO ex Su^SMM'^mm' (103) 

Therefore, collecting results we see that the operator O^^^ applied to the function 
l/(x%P*,i) in the closed case yields a term proportional to: 

■ 5ab5mn5{w -w')5w5rnm' (104) 

• '5AB(5AfA'(5(lt; ~ t«')'5n"^mm' (105) 

Where the constant of proportionality may depend on labels other than A and 
M and we have used the definition of K,q (j67| . The terms which depend on A and 
M are important as we shall seek to operate on ([50| with C"^^) which, via the 
directional derivative, links Qma to basis functions with different labels (for instance 
M + 1, A). Terms which depend upon the value of harmonic labels k,m,w however 
will be common to all terms and may be ignored if nonzero. 

For the K < case, the above argument may be repeated, relating derivatives of 
the harmonic functions to successive derivatives projected along (A )°. The operator 
O'^^^ applied to a function V in the open case yields a term proportional to: 



---(W)"<.i^^^ni^-"<--'(n<-^ 



l)i-ip- n+l)) 



2 

■ SabSaInS{w - w)5{p - p)5rnm' (106) 

= (2A+1)(A-M)! ^"^^("'"'^^ l^n "^J 

• SabSmnS{w - w')5{p - p)5jnyn' (107) 

For the K ~ Q case, the above argument may be repeated, relating derivatives of 
the harmonic functions to successive derivatives projected along (A )°. The operator 
C'^' applied to a function V in the open case yields a term proportional to: 

O^^'^V « .^^^^^±^^l^VMA{w,p,t) ■ 6abSmn6{w - w')S{p' - p)5^^, 

We now consider the effect of applying the polar projection operator, so as to 
pick out moments labelled N' and A' , onto the directional derivative d73l) i.e. 

O^P^P'ViV{x\P\t) (108) 
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A straightforward calculation shows that 

N' 2 2 

0^^'>P^\7,V{x\P\t) ex n«'EEQiv'^-i+«,A'-i+^4'^-i+o,A'-i+, 

q=0 a=0 13=0 

■"^N'-i+^A'-i+p ^2A' + 1)(A' - TV')! 

where again functions that arc common to both the results of ©'^'^(x^ _P% t) 
and C'(^^P*ViT^(x*, P*, t) are absorbed into the proportionality. 

4-. 2. The axial temperature perturbation 

We propose that for a given curvature K, the axial temperature perturbation field 
G'^"^''(x%P%i) may be decomposed as follows: 



e(^Hx\P\t)= J2 E 4/A(«^'fc)<]iKfc'^)<i (110) 



A.S)<M<A k,m,w 

where 

Awz r>A / ^\/^ab..M /^c „,de\T-7 C^ 

B^^\iz,ro,k,m,a,pn= ' ^^^^"^^ ^J^^ ^^"■■"^- (111) 

The difference between the axial and polar case is encoded in the appearance of 
the tensor tj"^-*^ instead of C"'''*'^. We shall call the former tensor an axial Chebyshev 
tensor. The properties of this tensor will be discussed below. As in the polar case, we 
shall seek to ultimately relate the moments Qj^^[^{w, k, t) to observables, choosing the 
functional form of the quantity c]y^jj(w, k) so as to simpliiy the appearance and 

4.2.1. Axial Chebyshev Tensors We define the set of axial Chebyshev tensors via the 
following scalar, vector, and recurrence relation 

g =0 (112) 

G"" - e,V (113) 

The tensors satisfy the following properties 

g.a.b.N ^ g.b.a.N (--^-^g^ 

Q-'^-^-^'lab = (116) 

g-^-^p, ^ g-^-i (117) 

g"''-^e,efc..e^ = sin(Ar cos-i(p^ee)) (118) 

where again Ca satisfies e°'e'^"iab — 1- 

An immediate implication of the identical recurrence relations obeyed by polar 
and axial Chebyshev tensors is that the axial equivalents of the functions Vna and B^jj^^ 
are affected by the directional derivative in the Boltzmann equation in an identical 
manner. Therefore we have that: 
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k.m,w A N=0 



Again we note that the axial metric source terms in the Boltzmann equation (|50p 
are simply expressed in the basis functions of the above expansion. As in the polar 
case, we now seek to find an operator that can isolate individual moments of a function 
which may be expressed in terms of a sum of axial moments. 

4-2.2. Axial projection operator Motivated by the similar formalism between polar 
and axial cases, consider an axial projection operator 0^^\ obtained by taking 
the operator O^^^ and replacing the polar Chebyshev tensor C'"''-*^ with an axial 
Chebyshev tensor g°-^-^' i.e. 

(119) 

where again O^'^) is simply related to O by swapping polar with axial Chebyshev 
tensors. 

The analysis of the effect of O'^^^ proceeds as in the O^^^ case. The significant 
difference is that one may no longer use the tensor W""^'^ ]^,j to related g°-^-^^ to 
p"'p''..p^^ . We instead introduce a new tensor Z"'''--^j ^^, which is a second variety of 
the operator 0, defined via the relation: 



nah-.M r7ah..M /icd..M /lonA 

By explicit calculation, we have 

ga = -{p+{A+)a-p-{A-)a) (121) 

I 

Then repeatedly applying the recurrence relation (|114l) . we then have that for 
M> 1: 

Gcd..M — -{2 ~ j5++..M(+)(A )cd..M (122) 

t 

-2^'-'p__„Mi-){A-)cd..M) (123) 

This implies that /3(^)+ — 2^^^^/i and (3(^a)+ — —2^'^^^/i, and therefore: 



■^° ■■ cd..M'^ab..M£p — -C^ "^ • V + + ..M^p(A )cd..M 

__ 2M-1 . V__..Mf,'(A(-'),rf..M) (124) 



However, it is readily seen upon evaluation that the integrals produce the same 
result i.e. a result proportional to 



Cosmological Perturbation Theory With Background Anisotropic Curvature 23 



'•'+'"' rWa(...M)fn.j^ 



{2A+1)(A-M)\ ^^' ' ' 'UJL «/ 

■ SabSmnS{w - w')6iii5jnm' (125) 

for ii' > 

i2A +l)iA- My. ^^'^^""^ "' ^> [^}^ '^"J 

■ SabSmn5{w - w')S{p - p')5mm' (126) 
for if < 

(A + uy. ^^^A), .. 

{2A+l){A-My^^'^^'^'^' ^ 

■ SabSmnS{w - w')5{p - p)5mm' (127) 

for K = Q. 

Similarly we have that 

AT' 2 2 

{A) (A) 



0(^)p^V.e(a;\ P\t) « n «' E E e)v'-i+.,A'-i 

q=0 a=0/3=0 



+[i'^N'-l+a,A'-l+p 



^(2-«,2-^) {A' + N'y 

■ ■''^N'-i+^A'-i+P ^2A' + 1){A' - N'y ^ ' 

Finally consider a function V'^^'^\x'^ ^P^,t) which is a sum of axial and polar 
components i.e. 

T/(-^^)(x^p^i) = ^FifMl+E^M>Mi (129) 

where the sum is over harmonic components and moments. 

It is readily checked that applying O'^^^Vj^jj^Bj^^j]^ = as it involves the following 
integral and sum: 

<^ J Y. rfV'c?'^(v++..Mf™'5(A(+)),d..M + V__..Mf;„'5(A(-)),d..M) • 

m' ,171 

(v++..Mf^(AW)='^-^ - V__..Aff^(A(-))='^-'^) (130) 

m' ,171 

= J J2 rfV'rf'^ (-v++..AfC'v— ..M^,^, + (-i)-('"+'"V__..MCv++..Mf;;'= 

in ' in 
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m' .m 

= (131) 

A similar argument implies that O^'^^Vpj^Bj^.j-^ = 0. Therefore the operators 
O^^^ and O'^^^ are suitable for picking out moments from functions composed of a 
combination of both. 

The Botzmann equation (|50p with no collision terms may then be cast as a system 
of coupled evolution equations for the moments Qj^J^{k,w,t) and Q}^j^{k,w,t) as 
follows: 

(i) Use the result ([7^ to express all metric source terms (and collision terms if 

appropriate) in equation (j50l) in terms of the functions Sj^^ and Bj^J^. 
(ii) Expand Q{x^,t) in terms of harmonic mode moments Qjyj^(k,w,t) and 
Q]^J^{k,w,t) using the expansions ([?T|) and (IllOp . 

(iii) Then use the operators O'^' and O^'^) to extract a particular term 
dQj^,\,{k,w,t)/dt or dQ)^,^,{k,w,t)/dt in ([50]) . This will be accompanied by 
terms in adjacent moments as dictated by the directional derivative of Q and 
(1721) as well as terms due to the metric. 



The presence of coUisional terms C j^ may be simply accounted for in the manner 
that the metric source terms are are if they are expressable in terms of the functions 

(P) (A) 

B\jj^ and B\,jj^ and do not violate the earlier assumption that the time component of 
the geodesic equation holds. 

5. Comparison of moment expansions in the FRW limit 

Having set up the Boltzmann equation, and found an appropriate projection operator 
to isolate polar and axial M and A components of the moment expansion, we now seek 
to motivate a particular choice for the functions c]^,j^. The polar case is of particular 
interest as it is likely that the source of the anisotropic, shear less background will only 
have polar perturbations (for instance the scalar field in Section |3]). 

Towards these ends, we consider the typical choice of moment expansion for scalar 
contributions to the temperature perturbation in the flat FRW case: 

e(A, <z, t) = ^(-i)^(2A + i)e^(g, t)VA{ii) (132) 

A 

where Va is a Legendre polynomial, /i is the cosine of the angle between the 
unit three dimensional co-moving wave-vector q of a plane wave perturbation (with 
magnitude q) and the photon 3-momentum. However, in the spirit of the 2-|-l 
decomposition of space, we may decompose the wavevector gq as follows: 



gq 



qPd, + qVl-P^k'^da (133) 

Where the 2-wave vector is unit with respect to the co-moving background 2- 
surface metric: jabk°'k'' — 1. This implies that 



f^ = p^q' = al3+ a/1 - /3 Vl - a2 cos(??) (134) 
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where rj is the angle between fc" and p° . 

Noting the group multiphcation property imphed by (matl), the addition theorem 
for matrix elements (addl) implies that: 

VA{[i)^ 2^ (-1) e '(:4^-p^^|M|(a)^|Af|(/3) 

M=l ^ '' 

^VA{a)VA{fi) 

+ 2 ^ CM(r/) J^-^P^(a)P^(/3) (135) 

M=l ^ + ' J- 

A temperature perturbation 0(Pi,x*,i) is expressed in terms of basis mode- 
moment functions PA{iJ)e.^^''^ . Using equation (|135p we have that: 

PA(M)e*^-^ = 7'A(«)PA(/3)e'^-^ 

+ 2 2^ ^ ^ PM(«)-PM(P)e (-fc) <- yabcA-ie 

M=l ^ '' 



n 

A 

I o \ ^ AnAncbi, \ ^ \ ) ' tjA f ^\ tjA / /d\Awz / i^\—M/-iahc..Al-s-i ck 

+ 22^^ e ^ 2^ TJ—J^^PM{Oi)PMKP)<i ("«) C. V abc.M^^ 

where we have used the Jacobi- Anger expansion: 

6*'^'''™'''^= ^i™J™(fcr)e™"* (136) 

m 

We see then that the mode-moment functions in the flat FRW case may indeed 
be related to the proposed mode functions S^.^j^ as follows: 



VAitAe^'^'' ^ Y^i''e"^*-VA{P)B^A^ 



n 

A 



(A-M)!^.„^^„^^ . , . iP) 



+ 2 E E Ji^^"^"'^^M(/3)eLl (137) 

The temperature perturbation Q{P\x^,t) around a spacetime with FRW 
symmetries admits the following decomposition: 

e(P^x^i) = y AEH)'^(2A+l)P^(P'<7,)e^(q,i)e''l'' (138) 
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Therefore using (|137p we may relate the moments 6a (Qj t) to the moments QnA: 
yielding: 



c'o/QoA = {-i)^i2A + l)eAi"e™^''7'A(/3) (139) 



and 



c^Pa^na = 2^^-^i-^)^{2A+l)eA{<l,t)^"^e^""^''P^{(3) (140) 
for N>1. 

(P) (P) f 

We make the following choice for the functions cjy^ = c)^^ in the flat case: 



(P)/_(-»)-^( 2A+l)(A -iV)! 
{A + N) 

Note that this choice in the flat case implies via (jl39p and (11401) that 






QoA =eA(q,i)«™e™*^7'A(/3) (142) 

Oata = 2eA(q, i)»'"e™*-Pj^(/3) (143) 

6. Free-streaming equation 

The linearity of the perturbed equations implies that we can isolate the impact of 
the anisotropically curved background itself on the evolution of inhomgeneities and 
isotropies of the photon distribution in the absence of perturbational source such 
as those due to the perturbed metric in the Boltzmann equation. In the absence 
of coUisional terms and metric source terms, each moment of 9(P*,a;*,i) obeys the 
following equation. 

= CnA&'nA + CN+UA-ie>N+l,A-l^"\N + 1, A - 1) 

+ cn+ua+iQn+i,a+iM"°{N + 1,A + 1) + cn^a^iQn,a-iM^^{N, A-1) 
+ CN,A+iQN,A+iM^°iN, A + l)+ cn^i,a-iQn-ua-iM^\N - 1, a - 1) 
+ cjv-i,A+iejv-i,A+iX'"(A^ -1,^ + 1) (144) 

The effect of the anisotropic curvature of space is via the appearance of the 
functions k^ in the matrix entries A4'^'^{a,l3) and A^*'^(a,/3). 

On the other hand, in the flat case K — 0, given an initial monopole perturbation 
Goo(fc, w;, io), the solution at a later time t is known: 

QNA{k,w,m,t) = 2JA{q{t-to))t"'e"^^''P^{(3) (145) 

where ja is a spherical Bessel function. This implies that (|145p is the solution of 

/ 

NA- 



(|144p when A' = given our above choice for the functional form of c 



Henceforth we adopt the following choice for cna'- 
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{-^)^{2A+l){A-Ny. 
^^^ "" iA , ATW (146) 

(A + Ny.KNKN-lKN-2--- 

- tJ^ (147) 

where Cjsfj^ simply means the choice for cna in the spatially flat case. It shall be 
seen that this choice simplifies the appearance of the photon temperature-temperature 
power spectrum (see equation (IB.26p ). Given this choice, we may write (|144p as : 

, . .{A-N)]^, k{A-N -2)1^ 

(^^ + ') (AT^®^^ ^ - "^+^^2 {A + N- 2)! ®^+^'-^-^ 

.k{A^N + 2)l^ 
"''~+^'2(A + iV + 2)!®~+^^-^+^ 
{A + N)\ ^ 

(A + l-Nl 
-^ iA-N)\ ^^'^+^ 

.fc (A -TV)! ^ 
+ '^"''2(A + iV-2)!®"'-^'^-^ 

k{A-N + 2y. , , 

+ ^^'2 (A + iV)! Q^-1-^+1 (148) 

where 8 at a is nonzero only for ^>0,0<A^<A. 

Care must be taken in the curved cases. When K < 0, the eigenvalue k'^ of the 
2-d Laplace-Beltrami operator need only take values greater than 1/4 to describe 
square integrable functions. This is not so in the flat case. Additionally, when 
K > 0, the function £^ is simply a constant when fc = and so, in defining the 
basis mode-moment functions Bna via gradients of this function, only iV = is 
permitted. In standard cosmology it is typically assumed that at recombination the 
photon distribution is described by a monopole and dipole Qq and 0i, each of which 
depend on the 3 dimensional g* only via its magnitude q. Inevitably then, contributions 
are considered from wavenumbers which lie entirely along the (9z)^ direction and so 
have no component in the orthogonal direction i.e. k = 0. We have seen above though 
that such a setup is complicated in the curved setup - if, for instance, 0oo depends on 
the combination q — ^JvP + k"^ then inevitably scales large compared to the curvature 
are probed as a constant value of q becomes more and more dominated by w; this 
does not appeal to an intuitive notion of the curvature being irrelevant at early times. 

By way of example of the effects of curvature, we initially restrict ourselves to 
the open case. We again restrict ourselves to an initial monopole perturbation, but 
now consider it to be sharply peaked at some value fcj ^ l-f'^l/4 with no z-dependence. 
Furthermore we solve for a perturbation n^r^ = ©JVA Ilq^o '^^/('*w'^a'-2'*^-4---)i 
and assume that A^ ^ 1 so that even though we are essentially restricted to 
contributions from wavenumbers far greater than those associated with the curvature 
scale, we may still have non-negligible deviations of kjv from unity. 
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{2A- 



.kN{A-N + 2y. 
''~ {A + N + -Z 
:kM jA-Ny. 
'' 2 (A + N-2y/ 



Y{A + N + 2y^''+'^^+' 



n 



N-l.A-l 






where we have introduced the variable fc^v = fc ^" " ^ " *'" . This equation for 

^7V-l^iV-3^iV-3'" 

the evolution of the 'moments' Una is, for each value of N, equivalent to the system 
in a flat universe and so admits the following solution: 

UNA{k, w, m, t) = 2JA{kN{t - io))i"e™*'=« P^(0) (149) 

and so 

e^^ = (^^4r-2^ir-4-) ^^.^(fe^(^ _ io))^"e^"^^«P^(0) (150) 

Now instead consider a case where instead the initial monopole perturbation has 
no dependence upon the coordinates on the co-moving 2-surface, now peaked at a 
value w = ki,. By the above arguments wc then have A'^ = and the Boltzmann 
equation takes the following form: 

{2A + 1)9^^ = wAQqa-1 - w{A + l)eoA+i 
There is again an analogue to the flat case, with solution 

QM^JA{w{t-t^))i'^P^{l) (151) 

Given that the initial monopole perturbation is sharply peaked at a value k■^, , 
we expect the Q moments in each case to be peaked at the lowest value of A for ja 
is maximized by its argument. For a function ja{x) this happens when A ^ x. 
We note then that this peak happens at A ^ k^,(t — to) in the latter case but 

2 2 2 

A ~ k^ ^A«A-2''A~4--- ^^ _ ^^-j jj^ ^Yie former case. 

Though far from exhaustive, these limited results seem to indicate that the notion 
of the position of an acoustic peak in A-space in the case of nonvanishing curvature can 
vary depending upon the value of N in question and the form of initial perturbations. 
This could have important observational implications. 
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7. The Boltzmann Equation for massive particles 

We now consider the Boltzmann equation for fields which can be described statistically 
in terms of a distribution function f'^ for massive particles. It's to be expected that 
this description is appropriate for baryons, cold dark matter, or massive neutrinos. 
The coUisional Boltzmann equation is then 

^ = C- (152) 

The time derivative df^ /dt is then decomposed as follows: 



df^ df^ df^ dx' df^ dp df^ da df^ dip'') 



dt dt dx^ dt dp dt da dt d{p<^) dt 



(153) 



As in the massless case, we expect dependence of /-^ upon x'^p"" and a as well as 
da/dt and dp"" /dt to be first order or great in smallness. Therefore up to first order 
in perturbations, the coUisional Boltzmann equation for massive particles is 

dJ^^dJ^dx^^dI^dp_^T ,,.,, 

dt ^ dx- dt ^ dp dt^^ ^^ ' 

We may define the field four-momentum in terms of the proper time r as follows 
P^ = dcc^/dr. We explicitly decompose the momentum as follows: 

P^'d^ = ^iP'dt+pP'd,) (155) 

P* EE aP* (156) 



P'd^ =ad,+pVl^^p^da (157) 

Therefore dx^ /dt — P^ /P*. The timelikc constraint P^P^ — —rri^, up to first 
order, is as follows: 



-m^= - (1 + y)(P*)2 + 2p(aE + {Ca + RaW^/l^^)P' 

+ p^ + U{1 - a^)p^ + FaV + 2(S'a + Ba)p"aVl-aV 
+ (l-a2)p2(^^^^2Q,b)pV (158) 

Where the metric perturbations here are the actual metric perturbations and not 
their harmonic modes. 

Up to linear order in perturbations we have: 



P* = y/p^ + m-^+p(aE + {Ca + Ra)pWl - "2) 
p^ .U . 2n F 



(77(1 - "') + 17"' + (Sa + BaWa^l - a^ 



\J p^ + rvt? 2 2 

+ i(l-a2)(X,b + 2Q„,)pV)-Vp^+^| (159) 
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As in the massless particle case, we restrict ourselves to cases where P* is not 
influenced by the coUisional terms that exist. From the time component of the geodesic 
equation we have then that: 

dP* P* d /P*\ ^^ 



Yielding: 



.. "d^ t)=-V^^'' (^«o) 



^ - -pH - ^ ^(aVa.P + 2e^aE' + e^ad.V + 2(1 - a')peU' 
+ 2pea^F' + a(l - a^)p^d^U + 2{e^ - p^)HaE)£ 



2e^C' - ApeaB' - eV - p^a^F)p''Va£ 



E ^^2e"' ^^ ~ "')p'^P"?5yV„V,V,£- 



VT-a 






- E "^^ "' (1 - a^y QP°Py VaVbVeg (161) 

where e^ = p^ + m^ . 

For a given species of particles, we can decompose /^(a;*, P*, t) as follows: 

f{x\P\t)^f^^,j^{p,t) + 5f{x\a,r,t) (162) 

where (5/^ is of first order in smallness. We propose the following harmonic and 
mode decomposition of the perturbation 5f^: 

5f{x\a,p%t)^ J2 E d^MAi^,k,t)4P^{w,k,t) 

k,m,w Afi<M<A 

(^ah..M\7 ck 

Jwz-pA(^\^ ^ ab..M<^„i 

■e t^J^.f(a)- 



{-kyt 

k,m,w Afi<M<A 

nah..M \-7 ck 

y y ab..MO^ 



.^^^zpA^^^y , "!m (163) 



where J"^^(zi;, k, t) and F\^j]^{w, k, t) are polar and axial moments of the massive 
particle distribution function. As such, the results from the massless case carry over, 
allowing for the additional terms that appear in (|16ip due to the presence of a a 
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nonzero mass. It is readily checked that the metric source terms appearing in (jl54p 
via ()16ip can be expressed in terms of the functions B}^J^ and B}fJ^. If this is also the 
case for the coUisional terms (having assumed already that they allow the use of the 
time component of the massive particle geodesic equation) then a set of time evolution 

(P) (A) 

equations for -F)^^ and -F^^l may be deduced precisely in the manner of the massless 
case i.e. by application of the appropriate polar or axial projection operator. 

The components of the energy-momentum tensor for the massive and massive 
fields may then be related to the distribution function and momentum 4- vectors P^ by 
standard methods, see for instance |26| . In principle, given a cosmological background, 
the collection of perturbed Boltzmann equations along with the perturbed Einstein 
equations are sufficient to calculate predictions for observables on scales where linear 
perturbation theory holds. 

8. Polar correlations in the photon distribution 

A very useful observable is the angular correlation function (Q^^^ {n^,^-.)Q''^\nl^-.) 
for perturbations to the photon temperature: 



e(^)(n|i))e(^)(n^2))) = EE («'''"'«'") ^™(n(i))^™ ("(2)) 



Im I'm' 



Where Yj^ are spherical harmonics and n^^ and nl^. are angles in the sky. 
Henceforth we will use the notation n = {9,(f)) to parameterize the angular vector 
n'. We now express the 'spherical harmonic space' correlation function {ai'm'OHm) ■ 



\^l' 7n' ^Im / 



^ fdn fdn' E E Y;Jin)Yi,in') 

•' •' N',A',N,Aw',w,k',k 

= ^ fdw fdn fdn' E E ^*" (^)^"' (^') • 



B^j^iiz = o,k,w,n){B\^lriz = o,k,w,n') ■ 

FJ,':\,^^ik,w,t) (164) 



where 



piP) (k ,,, f\ = ^(P)*iP) ^NA'^N'A' j^n N 

^N'A'NAy'^^'^^'') — ^NA'-N'A' Ig/f; y^)\2 'K'^t^) 

and we have defined the dark matter power spectrum {5(w' , fc', m')5{w, fc, m)) = 

V{k,w)5{w — w')Sk'kSmm', where S{k,w,t) is the dark matter overdensity and 6kk' is 

either an appropriate delta function or Kronecker delta symbol depending on whether 

( p) *fp) 
k is discrete or continuous. The step to express the expectation of 0)vA®JV'yi' ^^ 

terms of the dark matter power spectrum is a standard one in perturbation theory 

around FRW backgrounds and is also appropriate here; its utility is to isolate the 

contribution to due to the initial amplitude and phase of primordial perturbation but 
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that has no dependence upon P* . One expects this lack of momentum dependence to 
be the case for the dark matter power spectrum. Assuming adiabaticity of primordial 
perturbations, ratio Q/5 is not expected to depend upon the initial amplitude of 
perturbations. 

After lengthy calculation (detailed in [Appendix B[ ), one may carry out the 
integrations over ft and ft' in ()164p as well as the summations over N',A',N, and 
A yielding the following expressions for {ai'm'O'lm)'- 

K >0: 

{ai'm'aim) =^mm'y / du;^(2L + l) "' "'i; V{k,w)e{l,l',m) 



K <0: 



2 J ^' ' ' \5{k,w)? 

(165) 

QiP)Q*iP) 

{avm'ai^) ^5mm'2Tr\K\ / dw / dAAtanh(7rA) "'^ "'^ r{k,w)e{l,l',m) 

(166) 

q{P)q*(P) 

{ai'm'aim) ^6mrn'2TT / dw / kdk "^ "'^ V {k,w)e{l,l' ,m) 

(167) 
where the time dependence of perturbational quantities is implicit and recall that 



K = 0: 



k^ = kk* = KL(L + l)mtheciiseK >0,k^ ^ kk* = \K\{X^ + ^)) in the case isT < 0. 



9. Statistical anisotropy 

The CMB temperature anisotropy is said to be statistically isotropic if the following 
property holds: 



einl^Mnl^^)) - /(n^(i)n(2).) (168) 

i.e. that the two point angular correlation function depends only on the projection 
of the direction in the sky n(i) along the second direction n(2). Decomposing the 
angular functions, as before, in terms of spherical harmonics we have that: 



e(n|i))e(n{2))) = Y. {ai'm'aim) Y:!!inl^)Y:^inl,^) 

(169) 



,,l,ra'l' 



If (a/'m"^/m) — Ci^w^mm' ^^d (7/ is independent of tti, we have: 
;e(n|,))e(n^,))) = ^Gy^'(n|,))y^(n|2)) 



in,l 



O/ I 1 



47r 



/Ki)n(2).) (170) 
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In general {aiim'0,im) niay not in fact not be diagonal in I and I' or may depend on 
m, indicating then that the condition for statistical isotropy is not fulfilled. A natural 
origin for this in the event of background anisotropic curvature would be due to the 
presence of an additional direction in the problem: {dzY- In the spirit of our 2+1 

approach, we decompose the angles in the sky as follows n(i) = («(i), i/l ~ '^n-i'^n)): 

n(2) = (a(2), ^/l — Q^?2')"'?2))- ^^^ condition of statistical isotropy may then be written 
as follows: 

0(n(i))e(n|2))^ = / (a(i)a(2) + y/l - a^^^y^l - a^2)^a)"(2)a) 

We expect there to be a limited statistical isotropy in that the dependence of the 
correlation function upon h?^, and n?2) will only be via the projection of n^s along 
n?2')j there being no other preferred direction existing in the co-moving 2-surface. If 
we allow {ail rn' dim) to take the following form: 

{avm'aim) = 6„im'Cw{m) (171) 

then we have that 

e(n^i))e(n{2))) - J2 Cw{m)6^m'Y;:!:{nl,^)Yi{nl,^) 

m^l,m' ,1' 

ex Y. Gr(m)e™(^-^^)pi(a(2))/^^(a(i)) (172) 

where 0i and 02 are angles made between fi'^^-. and h%s with respect to a fixed 
direction in the co-moving 2-surface. Clearly though (|172p depends only upon the 
difference between these two angles and so only upon the projection of n'K-, along 
n?2')- We see then that the functional form of the correlation functions (|165p . (|166p . 
and (|167p in being diagonal in m and m' reflect the limited statistical isotropy of the 
problem. Clearly in general though, these functional forms need not be diagonal in 
I and V nor independent of m in the event that diagonality may exist. These are 
observational signatures of a preferred spatial direction in cosmology. 

One possibility is that this preferred spatial direction is encoded entirely in the 
primordial power spectrum i.e the subsequent evolution equations are not sensitive to 
a preferred direction. This may happen for instance if there is a period of primordial 
anisotropic cosmic inflation giving way to an isotropic FRW expansion afterwards. 
This possibility has been examined in detail in [7, and [8; . We may recover their results 
for the angular correlation function by using the result (jl67p along with equations 
(|142p and (I143P which relate the two-index moments Qma to the more familiar 9^ 
for scalar perturbations in the event that the background spacetime is spatially flat 
with FRW symmetry: 

(ap^'a;™) = 6rr.m'27r J dw J kdk^0^{q,to)r{k,w)P^iP)P^iP) 

I i2l + l){2l' + l)il-my.il'-my. 
V 167r2(; + m)!(r + m)! ^ ' 
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We now write the wave-vector space integral in terms of the magnitude of the 
3-dimensional wavevector q: q — \\/w^ + k^\ and /3, the projection of of g' along the 
z direction, yielding: 

{av^'Ohn) =^™m'27r fdqq^ f d(3^^{q,to)V{q,l3)Pi;^{l3)P^{p) 



\6\ 



'{2l + l){2l' + l){l-my.{l'-m)\ 
167r2(; + m)!(r + m)! 



(174) 



This is identical to the result of [5]. The ratio 0/0j*/|(5p is taken to be 
independent of the angular quantity /3 and so statistical anisotropy may occur solely 
from the following angular integral: 



Only if V{q,l3) = 'P(<?) is the equation ()175p independent of m and diagonal in I 
and /', therefore any angular dependence present will yield statistical anisotropy. The 
period anisotropic inflation considered in [7] and [8] corresponds to an example of a 
background where a ^ h and K = Q. It has been argued pi that one should expect 
the power spectrum V{q,/3) to respect the symmetry of the background spacetime. 
The absence of a dependence upon the unit wavevector fc" in the co-moving 2- 
surface is a reflection of this. Additionally though, there background has a reflection 
symmetry {z, x") — )• {—z, —a;") which suggests that one should additionally insist that 
P{q, /3) = Viq, — /3). This restriction implies that (|175p is only non- vanishing when the 
difference between I and I' is an even number ^. It is interesting to speculate whether 
this is expected to hold for the results (|165p and (|166p which arise from anisotropically 
curved backgrounds with the same reflection symmetry. 

10. The Inclusion Of Shear 

The preceding formalism has assumed that the cosmological background is shearless. 
There are advantages to these scenarios as they avoid constraints typically associated 
with shear such as those due to the anisotropy of the CMB present in the background 
and constraints on the anisotropy of the expansion rate. As discussed in Section [31 
shearless solutions with anisotropic curvature may exist only in the presence of a new 
matter field which allows the background solution to exist via the particular form of 
its anisotropic stress. In the absence of such a matter field, models with background 
anisotropic curvature concomitantly possess shear ^, and therefore a more general 
analysis of the effects of background anisotropic curvature should allow for background 
shear. Recall that allowing for shear, the background metric can be written by as: 

ds^ = -a^ (dt^ + —dz^ + -fcddx^'dxA (176) 

In the presence of shear (i.e. b{t) ^ o,{t)), the co-moving metric is clearly no longer 
static. However, with this choice of co-moving metric we may as before decompose 
metric perturbations in terms of eigenfunctions, eigenvectors, and eigentensors of the 
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the Laplace-Beltrami operator associated with the static metric jcd- For instance, a 
function may be again decomposed as follows: 

y(i,z,a:^)=^y,(i,z)f^(a;'^) (177) 

l,ni 

It is then desirable, as before, to decompose the functions yk(t,z) into harmonic 
components and then obtain field equations for modes which depend only on 
wavenumber labels and time. 

A choice for such harmonics are the eigenfunctions 6^,(2:, t) of the Laplace-Beltrami 
operator on the co-moving 3-space projected along the z-direction, which are thus 
defined by the following equation: 

-^jd^_e^iz,t) = -wh^{z,t) (178) 

Where the w'^ are constants. Solutions to this equation are then simply plane 
waves but with time-dependent frequency w = {b/a)w, i.e. 

ykit,z) ^ ^ J Mt)ykAt)e'''^'^' (179) 

Therefore one may decompose a perturbed tensor such as the Einstein equation 
in terms of generalized sums of the functions e„, multiplied by the eigenfunctions, 
eigenvectors, and eigentensors of the the Laplace-Beltrami operator associated with 
the static metric 'jcd- Care must be taken with partial time derivatives of (jl79p as the 
integration measure and 'wavenumber' now have a time dependence. Consequently, 
the components of the perturbed Einstein tensor 6G'\ will involve additional terms 
due to the presence of shear in the background and its effect on the definition of 
harmonic components. 

If the shear is sufficiently small, it may be regarded as a small homogeneous 
perturbation to a shearless background spacetime. For instance, in the closed case 
this may be realized by considering the polar perturbation F = F(t)5{w)dLo, for 
which the collisionless photon Boltzmann equation becomes: 

(eo2 + ^F)' = (180) 

(eoo + ^^)' = (181) 

The contribution to 8ooj assumed to have no spatial dependence, can be 
absorbed into the definition of the background temperature, leaving a time-dependent 
contribution to {aimdi'm') = Ci{m,t)SwSi2Smm'SmO, i-e. a quadrupole contribution 
to the CMB. Note however that the perturbations to the CMB are not statistically 
isotropic as the function Ci depends upon m. 

In treating the presence of shear as a background quantity, further complications 
emerge. Consider decomposing the momentum of a photon propagating in the 
background spacetime as follows: 

We decompose the particle momentum P^ as follows: 
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p^'d,, = p'd, 



1'^* 



a 



k + 1{{1) pad.+pVl-a^p'dc) (182) 

P'dt + P'dr) (183) 



The null-norm constraint immediately enforces P* = p if it is required that 
IcdP'^p'^ = 1. The geodesic equations then imply the following relations: 



da 
dt 


= a{a^ - 


-1) 


fb' 


a' 
a 


dp" 
dt 


= 









it — n"5+(^-")7j (^^4) 

(185) 

(186) 

Therefore one can no longer assume that ^ = at the background level, and this 
must be accounted for in the Boltzmann equation along with the background angular 
dependence of / implied by the presence of shear. However, it seems that the expansion 

(P) (A) 

of temperature perturbations in terms of Bj^j^ and i3)vA ^'^^^ remain appropriate, and 
that the expressions for the angular correlation function ()165p . (|166p . and (jl67p will 
still apply- the difference from the shearless case appearing in a generically different 
dependence of the moments Qna upon the labels w, k, and t. 

11. Conclusions 

In this paper we have presented a formalism for the analysis of cosmological 
perturbations to an anisotropically curved but shearless background. In particular, 
a harmonic decomposition of spacetime dependent fields has been introduced, as 
well as a moment decomposition of the dependence these harmonics may in turn 
have on the direction of particle momentum. One then obtains via the massless 
and massive particle Boltzmann equations a system of coupled ordinary differential 
equations in conformal time. These equations are coupled to one another by the 
presence of adjacent moments of an identical species, possible collisional terms, and 
to the spacetime via metric source terms. In turn, the Einstein equations describing 
the evolution of metric perturbations are coupled to the matter fields via their stress 
energy tensors which can be determined from the fields' collected harmonic moments. 
Given a set of initial data, these equations are sufficient to describe the evolution of 
cosmological perturbations. 

Of particular interest is the temperature-temperature power spectrum. Explicit 
expressions for the angular correlation function for polar perturbations have been 
determined in Section [51 and the conditions under which they represent statistical 
anisotropy have been discussed. In particular, these expressions have been written 
in a manner which seeks to separate statistical anisotropy which is present via 
initial conditions (primordial statistical anisotropy) from that which ensues from the 
evolution of perturbations on an anisotropically curved background and in the presence 
of perturbations to which fields produce this curvature (which may be termed emergent 
statistical anisotropy). Indeed, from our formalism we have recovered antecedent 
results regarding statistical anisotropy of temperature fluctuations induced by an 



Cosmological Perturbation Theory With Background Anisotropic Curvature 37 

anisotropic primordial power spectrum in an FRW universe (an example of primordial 
anisotropy, see [Z],IH]) as well as statistical anisotropy induced by a background with 
anisotropic curvature and anisotropic expansion rate (emergent anisotropy, see [S]). 

Although in principle the formalism may be used to determine the evolution 
of cosmological perturbations, there are yet several steps to take in order to be 
able to accurately assess what kind of observations to expect in a universe with 
background anisotropic curvature supported by a field such as that discussed in Section 
[3l Firstly, it is important to provide a careful treatment of primordial power spectra 
in the presence of anisotropic curvature, where intuitively it would seem that the 
effects of curvature should be unimportant for the co-moving scales probed today. 
Secondly, the results in this paper have been presented in a general gauge and the 
existence of particular gauge choices which may offer considerable simplification has 
not been determined. Thirdly, polarization in the photon field has not yet been allowed 
for in calculations, and possible collisional terms (for instance Thompson scattering 
between electrons and photons) have not yet been explicitly included in the Boltzmann 
equations. Additionally, it would be desirable to extend the formalism to the case 
where shear is present in the background. The modifications to the formalism that 
this should involve have been discussed in Section [101 
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Appendix A. Special functions and addition theorems 

In the decomposition of spacetime dependent perturbational functions, vectors, and 
tensors we have made use of the eigenfunctions £^ of the Laplace-Beltrami operator 
V^ on the co-moving two surface. We now discuss properties of these functions in 
each curvature case in more detail. For simplicity we discuss eigenfunctions of the 
Laplace-Beltrami operator V^ which is associated with the co-moving surface with 
metric 'jcd (see Section [1] for discussion of this quantity). 

Appendix A.l. Positive curvature case 

In the positively curved case eigenfunctions 5^(5:°) of the Laplace-Beltrami operator 
V^ are the familiar spherical harmonics F^, which satisfy: 

v2y,U^» = -i{i + i)YUo.^) (A.l) 

-.^^^^y:(M) (A.2) 

where / is an integer taking values between 0, 1, 2...oo and ~l < m < I. 

Appendix A. 1.1. Connection to representations of the group SU(2) Recall that the 
group SU{2) has a realization as the collection of all 2x2 matrices Mab with complex 
coefficients such that Mab = ^^^ba ^^^^ dct(A/) = 1. It may be shown [27 that a 
generic element g of the group may be parameterized by three numbers 4>, 9, ^p such 
that M{g)AB takes the form: 



cos(|)e*(^+'^)/2 zsin(|)e*(^-'^)/2 
^ sin (I) e'(^-^)/2 cos (|) e-*(^+^)/2 
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where < < 27r, Q < 9 < 27r, ~2tt < i/j < 2Tr and the identity element 
e of the group is the matrix MjiB{9,(j>,4' = 0) = Sj^g. With this realization, the 
group multiplication operation which takes two group elements gi and ^2 to produce 
another element gig2 is simple the process of matrix multiplication: MAB{gig2) = 
MAcigi)^^^ ^103(92) and so, for instance "012 may be fovmd in terms of ^i and tp2 
etc. 

We now adopt the following definitions of |27) : 

• Let V be a Hilbert space. A representation of a group G with representation 
space V is a homomorphism T : g — > T{g) of G into the space of bounded linear 
operators on V. From this it follows that: 

T(5i)T(52)=T(5i52) (A.3) 

T(5)-i =T(5-i) (A.4) 

T(e) =1 91,92,9 &G (A.5) 

• The representation T is reducible if there is a proper subspace W of V which is 
invariant under T. Otherwise T is irreducible. 

• Let {vn} be an orthogonal basis of V with inner product («,„,««) = Smn- The 
matrix element T,„„(g) may then be defined as follows: Tjnn{9) = ("^m, T((7)wm). 
If T^^{g) = Tjnn{g~^) then the representation is said to be unitary. 

It may be shown ([^) that for the group SU(2)), there're a set of irreducible 
representations T labelled by the number I which is an integer ranging taking values 
0, 1, 2...C)o, with basis vectors e*™" {—I < m < 1,0 < a < 27r) and accompanying inner 
product (e™",e™") = ^ /e-*™"e*""da. In this case the matrix elements TJnl take 
the following form: 



T^li9) = i"-"e-'(™^+"'^)pW^(cos e) (A.6) 



where 



T?(i) 



jl - ny.jl + m)\ 
{I + n)l{l - my. 



{I - r])'"^ {1 + Tj) 
2™(rn-n)! 



1 + r. 



X2Fi(/ + m + l,-/ + TO;TO-n + 1;— -^) (A. 7) 

where 2^1 is Gauss's hypergeometric function. The above expression applies when 
m — n > 0. If m — n < then one must replace m and n by —m and —n respectively. 
This function has the following symmetry properties: 

-PiLiM = (-i)"+"nl2.(^) (A.8) 

^^l(^) = i-iT'-'-P^-l.-niv) (A.9) 

T'^K^) = P^-l-nAri) (A. 10) 

It follows immediately that 

P^o)(cos0) = Py)(cos0) (A.ll) 
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where Pi{cos9) is the famihar Legendre polynomial. Furthermore we adopt the 
following definition of the associated Legendre function Pljt : 



Pa)iv) = i-i) 



, T{l + m+l) 
ni\T{l — m + 1) 



1 — 77 



1 + 77 



V2 



X2Fi(l + /,-/;m+f;iy^) 



for TO > and 



no(^)-(-ir r!!!I^tn ^(o"(-/) 



r{l-m 



if m < 0. 

Therefore it follows that 



Pmiv) 



{n + iy. 



(/ - n)l 



1/2 



,(0 

-n,0 



where n > 0, and hence: 



Ai) 



n;o(.9) = *-" 



4tt 
21 + 1 



YnAg) 



where Y"^{g) is the spherical harmonic function, defined as follows: 



(A.12) 



(A.13) 



(A.14) 



(A.15) 



' ^^ y 47r(/ + m)! ' ^ ^ 



(A.16) 



The group multiplication property (jA.3|) produces the following result for the 
matrix elements: 



Tm9i92) = J2 T^'\9i)mp6P'T^l^ig,) 



(A.17) 



As the representation is unitary, we have T^h (gi) — Tj^a (Si ) ^^'^ ^"^ ^^ have: 



q=-l 

Applying the definition (|A.15[) we then have that: 



(A.18) 



47r 



21 



- j yo'((^» = 0) = Y.Pol'^aiWT^'ligi) 



(A.19) 



E 



All 
21 + 1 



y::{g)yL{9) (A.20) 
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This result is known as Unsold's theorem and will be useful in our calculations. 
An additional summation that appears is the following: 

Y^ i(-)-''-»v,b..„yAii(+)''^-"v„fc..„r;/' (A.21) 

We now show that when A^-^^"' is defined in using eg and e^ as basis vectors, 
the summation (jA.2ip can be shown to be an example of the addition theorem. The 
matrices Vmn satisfy the following recurrence relations: 

arj y'l — r]'' 

VT^ivJil - ^^^J=4n'^l = V(l + n){l-n + l)vli^,„_, 



We write these equations in terms of matrix elements: 

(^^^^^ "7T^ ^ 'MJ ^™" " -V{l + n)il-n+l)ze-^^T!^l_, 

We now consider the tensor Xbcd..m = '^bcd..mf where / is an appropriate 
function. Now consider taking the gradient of this tensor and projecting all indices 
along the vector field i'^^" = \K\-^/^A^-'>'' = {{ee)" + i{e4,)''). Therefore we have 
the identity: 

- X,,..,„i(-)'^V,i(-)''^-™ (A.22) 

The gradient of A^'''" is: 

Vbi(-)'^ - ^{dM^iieeT + (e^T) (A.23) 

It follows then that V^A^-^'^^-™ is: 

y^j^(-)bc..m ^ ^i--("-i)Vai(-)j' (A.24) 

j 

= ie^),^^J2 A^-^--^"^-'H^^{eoy + (e,)-'KA.25) 
V2smt^ 

= " '(^^^-^iS ^ A^~^--^--'^A^-'>' (A.26) 

j 

= -zm(e^),^i(^''-'" (A.27) 

sm 
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where the last hne fohows from the symraetry of A'"^--"^ with respect to the 
interchange of all indices. Therefore: 

"^C0S6> 2(-)bc.,mv 
7^~. 7. ^bc.m 

V2sm9 
similarly it may be shown that 

A(+)-A^+')b--^VaXtc..,n = A^+^-Va{A^+'>''-"'Xtc..,n) 



V2sm9 

l{ + )ay TOC0S6> \ ^(+)(,c..niv 

V2sm9j 
Additionally, on some function y we have 

d+y = ^{dg + -^d^)y = A(-^-Vay (A.28) 

V2 smt^ 

d^y=^{dg-^d^)y = A^+^'^Vay (A.29) 

Therefore, projection of the gradient of Xbcd..p along the vector ^(+)° for all 
indices is equal to the differential operator I'+(p) = {1/V^){d0 + {i/ sin 9)84, — 
p{cos 9) / {sin 9)) applied to the projection of Xbcd..p along the vector ^4^+)° for all 
indices. 

However, recalling that 77 ~ cos 9, we similarly see that the recursion relations 
imply that the matrix element index n is raised by operating on Tmn with a differential 
operator proportional to I?+(p): 

V+{n)T^l = ^^il-n){l + n+lWT^l„^, (A.30) 

Now recaU that Y„j{g) = ^"^^(2? + l)/47rTl'^o(.g). 
This implies that 



A^-^'^'-PWabcpYnl^ l[V+{q)Yj^ (A.31) 

= []p+(g)z"J_^ri„o (A.32) 

9=0 



V2y V 4^ -"'p 
p 

n ^J{l + q){l-q+l) (A.33) 

g=0 
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Note the close relation of these projected derivatives to the spin- weight spherical 
harmonics [5S],[21]. Then rearranging (IA.33P we have that: 



(^) 



, -n—p—ip<S/ I JTT 

' ^ ^' 21+1 ~, 



T^') =W£Z V i!li 4(-)a^-lplv , . .yi (A 34) 

-"■p ttIpI rrr, — \n r~rT ^abc..\p\^n v^.oi; 

Ugio vC + ^Ji^-^ + i) 

The addition theorem for matrix elements immediately implies then that: 
47r P'P^ ^ 47r 

-2p 



(^)" 



m— ^ 

-^^ V ab-.p-'- rri-'^ ^ ab..pJ- ni 



Appendix A. 2. Negative curvature case 

As in the positively curved case, we characterize perturbations in terms of 
eigenfunctions of the Laplace-Beltrami operator V^ now defined on the upper sheet 
of the two dimensional hyperboloid of two-sheets possessing metric -jab- We shall call 
these functions Z^, defined as solutions to the following equations: 

y^ziie, 4>) = - ~ezi{e, 0) = - (^i + a^^ zi{e, 4>) (a.35) 

■dZj^{9,<j)) ^rrk (Q i\ /A na\ 

'd4> " "^^m(^' W {AM) 

Writing k'^ = k ■ k* , solutions appear in three branches. The first covers values 
k^ > 1/4, for which k takes the value k — —1/2 -I- iX where < A < oo. The second 
covers values k = —1/2 + i\ where Oi < A < — (l/2)i. For these two cases the label 
m satisfies m S Z. The third branch covers the case k = 0, allows only m = and is 
given by Zq — constant. 

Note that a square-integrable function f{9) on the interval [0, oo) can be expanded 
as follows: 



f{0)^ c(A)Qoo'^*^(^)^tanh(7rA)dA (A.37) 

Jo 

where 

c(A) = / f{e)Q^fj^^'^{e)smhede (a.38) 

Jo 

i.e. to describe square integrable functions on the hyperboloid we need only 
consider the first branch < A < (X), up to a constant component to the function. 

We now outline how eigenfunctions of the Laplace-Beltrami operator on the upper 
sheet of the hyperboloid of two sheets are related to matrix elements of a unitary, 
irreducible representation of the group SU(1,1). This group is the collection of all 2x2 
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matrices Nj^g which are unitary and have determinant equal to unity. An element g 
of the group may e parameterized as follows: 



cosh (I) e'('^+^)/2 sinh (|) e'i^-^^^ 
sinh (I j e^('A-«/2 cosh (|) e-*(^+^)/2 



where < (p < 2tt, —2tt < i/j < 2tt, < ^ < oo. By inspection the identity element is 
associated with (p — tj; — t — 0. 

It may be shown gS] that there exist irreducible representations T tj^ e) where 

where k = iX — 1/2 (0 < A < oo) and e may take the value or 1/2. The basis for each 
of these representations may be taken to be eigenvectors e""" (— oo < m < oo,m e Z, 
< a < 27r) and accompanying inner product <^e*^", e*'*") = ^ J e^'^^e^^^da. As 
before, matrix elements T^j, may be defined as follows: 

TM(g)e-^io^ ^ T^^\g)e-'''' (A.39) 

b—-—oo 

they are found to be [25] : 

T!,i\9) = e-('^+^)<^-'(''+^)'/'Q^+,,,+,(cosh(0) (A.40) 

For a > b we have: 

r ^ r(fc-6+l)(cosh(|))2fc(tanh(i))"-^ 

''''^^°' ' {a-b)\T{k-a + l) 

X2Fi{-k-b,-k + a;a-b + 1, tanh^(-)) (A.41) 



For a < 6 it is necessary to replace a and b by —a and —6 respectively. The Q^j, 
satisfy the following symmetry relations: 



Q:;„(coshO = Q*:™,-„(coshO (A.42) 

QL -Tl^^^^^riSl^^SijcoshO (A.43) 

r(fc + n + l)r(fc-n+ 1) 

The addition theorem for the matrix elements is then: 



+ 00 

T^b\9i92)= J2 Tii\9i)5'''T^l\92) (A.44) 



C— — OO 



have: 



Agam, as the representation is unitary, we have T^\gi) — T^^ idi ^) ^^*^ ^^ "^^ 



+00 

7^i?^(</' = ^ = e = 0)- Y. Ti^H9iW''T;^''\9i) (A.45) 
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From hereon in wc shall look exclusively at the case e = 0. It may be checked 

; the function T^^^^ ^ is an eigenfunction c 
upper sheet of the hyperboloid of two sheets: 



that the function T^^^^ ^ is an eigenfunction of the Laplace-Beltrami operator on the 



VX:-Hcoshi9)) = -[^ + xA r:VMcosh(0)) (A.46) 



where 



~2 _ _^ cosh 6* d 1 d^ 

The matrix elements T_^q = e^"^'^Q^{coshd) = Z^ may then form a basis for 
a harmonic expansion. The addition theorem for matrix elements then implies that: 



Zl{{0, 0) = 0) = 1 = ^ Z*^{cosh{e))zi{cos\i{e)) (A.48) 

— oo 

This is an analogue of Unsold's theorem theorem for the case of functions on the 
hyperboloid of two sheets. Another quantity that appears in the calculations is the 
following summation: 



m— — oo 



i(-)"'-^V„b..pZ£i(+)"''-PV„b..pZi (A.49) 



Analogously to the positively curved case, we will seek to calculate this quantity 
by relating projections of V ab..pZ'^ along the vector A^~^°- . 

In the hyperboloid case, expressions for A*-*'" in terms of unit basis vectors eg 
and e^ are identical, however now (e^)" — {1 / smh9){d(f,Y . Explicit calculation shows 
that: 

V.i(-'" = ^°^^^g (e0)b(-»(e.)° + (e,)") (A.50) 

Following the calculation of the positively curved case, we find that for the tensor 

■^hc.m = ^ bcmj 

m cosh 6* T(_)bc..m^ 



\/2sinh0 

V V2sinh^; 

similarly we have that 
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m cosh 61 ^(+)fcc..my 



ySsinhfl 



Additionally, on some function y we have 

9+y = -^i^o + -^d,)y = A^-^^VaV (A.51) 

5-.= ^(a.--^a,), = iW"V.. (A.52) 

Therefore we have the recursion relations: 



A^-^'''-"'+'Vab...m+i = I?+(m)i(-)°^-"V„fc...™ (A.53) 

i(+)"'-"+'Vafc...™+i = 2?-(m)i(+)''^-"Vafc...™ (A.54) 
where 

. X 1 / r, * „ cosh \ , . s 



It follows then that 

i(-)'^'-^V,fa..pZ£ = f[v+ip)zi (A.57) 

q=0 

Prom [3S] we can read off the following recursion relations: 

d ^«-^-«^^^^2^. = (^-^)QU (A.58) 

Qab^ik + b)Qi,_, (A.59) 



de V sinh 6 
d f h cosh 9 — a 



je V sinh 61 
or in terms of the matrix elements 



therefore, recalling the definition of T^l in terms of constituent functions 

we have that 
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d_ 
09 



— b cosh 6 
sinh 9 



ri = V2V+T 



+ -^ab 



(jA-1/2) 



= V2(fc - 6)e^*rf,,+i 



and therefore 



lA^'A 1 



^.= bj niT^^-<'>^'' 



q=0 



- -i n^^.(^)^'^. 



V2 
V2 



q=f)k + 1 



'ip^i 



\P\ 

n 



=0^ + 9 



i^-^'"-l''IVfc,..|p|Z^, 



(A.62) 
(A.63) 

(A.64) 
(A.65) 
(A.66) 



Then applying this equivalence to the addition relation for matrix elements, we 
have that: 



, 1 \ -P 

rp{iX-l/2) _ ' ^ 



n 



g=0 ^ + 9 



+ CXD 



1 



J2 S'fA^-'>'^-PVbc..pZ'lfA^+'>''-PVbc..pZ*_\ 

(A.67) 



e— — oo 



Appendix A. 3. Flat Case 



Finally we consider the flat case. The appropriate group here is the group IU(1), 
elements of which can be represented by the following matrix: 

e*" z 
1 

where < a < 27r, z G C The basis for each of these representations may be 
taken to be eigenvectors e*"*" (— oo < m < cx),m e Z, < a < 27r) and accompanying 
inner product <^e*^",e'^") = -^ J e^'^^°'e^'^°'da. As before, matrix elements T^^ may 
be defined as follows: 



+00 



they are found to be [! 



b— — oo 



Ti'^^(ff) = i"-"e-"^-'("-'")'^J„_™(fcr) 



(A.68) 



(A.69) 



where Jn-m{kr) are Bessel functions of the first kind. We define the following 
functions: 



J^ ^ T^_„ = z-"e»* J_„(ftr) 



(A.70) 
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These functions are indeed eigenfunctions of the Laplace-Behranii operator on 
the Euclidean 2-planc, and satisfy 

V2j'W(fcr) =-fc2^«(fcr) (A.71) 

-^^^^^ = nji'Hkr) (A.72) 

The Bessel functions satisfy the following recursion relations: 

^ + ^ j Mz) - J._i(z) (A.73) 

^ - -) JA^) = - J,^+ii^) (A.74) 

dz z J 

or in terms of matrix elements 



d_ {n-m) \ ^(,fc) _ ,i.^.VT(^fc) 
d9 






in I " ^ ^m-l,n (A. 76) 

Again it is useful to define vectors A^^^"- = {l/^/2){{egY T i{^<tiY) where 
eg = dg and e^ = {l/9)d^ i.e. they are unit vectors with respect to the metric 
ds^ = dO"^ + 6^d(j}^. Explicitly we have that 

V,i(-)" = -l-(e^),(-z(e,)" + (e^)") (A.77) 

Following the calculation of the positively curved case, we find that for the tensor 



^bc.m 



™ i(-)b-mx. 



V20 



A(-)-Va -^] i(- 



)bc..m 



V20J 

similarly we have that 



Xbc 



Ai+^-A^ + ^bc-rn^^Xbc.rn = A^+'^^W a{A^ + ^'''-"' Xbc.rn) 



^ A(+)bc..rr.x. 



V20 



Additionally, on some function y we have 
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d+y = ^{dg+'-d^)y = A^-^''Vay (A.78) 

a_?y = -^(ae-^a^)?7 = i(+)"v,iy (A.79) 

Therefore we have the recursion relations: 



i^-^"''-'"+'V,b...„+i = P+(m)i(-)"^-"V„fc...„ (A.80) 

i(+)"'-"+'V„fc...„+i = P_(m)i(+)"^-"V„fc...„ (A.81) 



where 



V+M^^lde+'-d^-"^) (A.82) 

1 

V2 



V^{m)^^(do-'-:d,-'-i\ (A.83) 



It follows then that 



i^-'""-^ V„fc..p£^ = n ^+ (P)^™ (A.84) 

g=0 



Therefore we have that 



rif-(^)'(^) ^^-"e-^^'^i'"-I^IV6,..HJ'!'j (A.85) 

Then applying this equivalence to the addition relation for matrix elements, we 
have that: 

/ , \ -p / , \ p +00 



2/ Vfc2 

^ ^ ^ ' e— — 00 

= 1 (A.8( 

In the case p = one recovers the familiar addition theorem for Bessel functions 
of the first kind. 

Appendix B. The calculation of polar correlations 



= JdnJdn' Y. E y^'mYi^m 

N' .A' .N.A w' ,w.k' ,k,p.p' 

B\^i{x\ k, w, n){B\^lr {x\ fc', w', n')Flj:\,^^{k, k', w, w', t) 

= ^ fdw fdn fdn' Y. E YrmYi'^' (^') 
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b\^1(z - o,k,w,n){B'j^ly{z = o,k,w,ii') ■ 

FJ,%j,Ak.w,t) (B.l) 

fn,^^(fc,z«,t) ^ c^Nl4'I'^f£^nk,w) (B.2) 



■^k,w 



where recall that here k and p are labels on the harmonic functions £* and ^^ 
and ^j, denote the appropriate sums/integrals over labels k and w. 

Iff we work in the preferred frame, then a particularly convenient choice is 
dCl — —dadx where x is an angular coordinate in the tangent plane at s some point 
in the co-moving 2 surface. 

We now consider the integral over x- The integrand will contain terms e"*™-^ from 
the spherical harmonic F,*' as well as terms C"'''"^ = M^"'' "^cd mP'^'^"'^^ contained 
within the mode functions Bna- It will be helpful to express e*™-*^ in terms of tensors 
formed from p'^ . We choose a local orthonormal basis (ei,e2) in the co-moving 2- 
surface, such that 



p'' = cosx(ei)'= + sinx(e2)" (B.3) 

= \{e'^ + e-'^'KeiT + ^(e^'^ - e-^'^)(e2)^ (B.4) 

therefore we have that 



e^x =p^((ei), + i(e2)c) (B.5) 

e-'^=pmei),-i(e2)c) (B.6) 

We claim that without loss of generality we can choose the coordinate x such that 

((ei)c + *(e2)c) = V2(A-), (B.7) 

((ei)c-*(e2)c) = \/2(A+)e (B.8) 

and so for a quantity e^*™-^ we have: 

g.lmix ^2^(A-)„,..|,„|p'^^-l™l (B.9) 

e-'^"'^'' - 2^(A+),b..|„|r''-l'"l (B.IO) 

The integral over x is as follows: 

For TO > we have: 

■rf^e-™'^C'^^-l^l=2^(A+),/..|„|T4^''^-l^l,..l^iy'dXP^/^'^-l'"l+l^l 
whilst for ?Ti < we have: 
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The above integral may be evaluated using standard methods, yielding: 

'^ -ab..|m| + |JV| . ^ (|m| + |A^|)! OTr^{ab..{\m\ + \N\)) 

P "^ 2(IH+|Af|)(((|m| + |M|)/2)!)2 ^ 







which is nonvanishing only for {\m\ + \N\) being an even number. 
In general then we will be left with terms like: 

(cdef..{\7n\ + \N\))Txr<^b..\N\ (A^),, iV^i^n^*^ ^12^ 

— 7 \^ )ef..\m\\'^ ^ ++..\N\^p\-^ )cd..\N\ 

+ 2l^l-iV__..|jv|fp'(A-U.|^,|) (B.13) 

Again we note that by the definition of the tensor W"' " \,^ lArh i* must be 

traceless with respect to contraction with ^'^^ across any two lower indices. Therefore, 
(|B.12[) is only nonvanishing if '^(''''■■(l^l+I^D) is contracted with one index of each 
component ')'^^ acting on one index of the W tensor and one index of the A tensor. 
This implies that \m\ = \N\. By a counting argument, there're 2'^Q\{Q)\ permutations 
of 7"^-2Q that fulfil this condition. Therefore: 

In total then we have that: 

V<,,..|jv|£,^y'dxe-*"^C"''-I^l = 2^7rV^^..|„|£,*<5|„||^| (B.15) 

Where — and + correspond to the cases m > and m < respectively. The case 
m = is trivially evaluated, yielding: 

fp^y'c'^''-l^ldx = 2^fp%l^l (B.16) 

We may immediately then evaluate the similar integral (again assuming m ^ 0): 

^ab..\N'\S;'' j dxe™''^C'^''-l^'l = 2^7rV±±..|™,|£:f <5|™,||jv'| (B.17) 
We next consider the sum over p: 



p 



2 2 7rV±±..|,m'|fp '5|„j/||7V'|2 2 7rVzpzp..|„j|£p(5|TO||7v| (B.18) 
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We first make the physical assumption that m — m' (and hence N — N'). 
Following the arguments in Section [HI this is a statement about the statistical isotropy 
of correlations in the co-moving 2-surface. This assumption simplifies (JB.ISP to: 

J2 2l"l7r2v±±..|„|£:f V^^..|™|fp^ (B.19) 

p 
We now may make use of the results of |Appendix A[ relating the derivative terms 
in (jB.lSp to matrix elements for group representations, where it was found that: 

2fc+ 1 , 2fc+ 1 



47r i"'ii"'i^ ' 47r 



r|m| 



/ . v++..|„|ypV — ..\m\yp 



It follows then that in the case K > Q: 

Y: 2l-I^^V±±..|™|£;'=V^^..|„|£,^ = ^'^ + ^^"'''''"' ^f[{l + q){l-q+l) 

p——k q—0 

Again we may directly use the results of [Appendix A| in the manner above, to 
find: 

+00 \m\ 

^ 2i™i^2^±±..i^|£:;'=v^^..i™|£:p^ = ^2|i^|l™l Y[{k + q){k* + q) 

P— — CO q—Q 

Additionally we must integrate over the a and a'. In the preferred frame there're 
two functions with dependence upon a: the associated Legendre function P^{a) 
appearing in the spherical harmonic Y*j{a,ip) and the associated Legendre function 
-^jv=lm|('^) appearing in iJjvA- ^^ identical combination appears for the dependence 
upon a' . The form of (jB.ll) dictates that we must integrate over the product of these 
functions, which is a standard integral: 

;' A' I 2(l' + m)\ 



for ?7i > and 



jpLP(L\da = (-I)l"-I jj^|^|j; /Pi<|da (B.22) 



TI<-)-(f^)^^' 



(B.23) 



[21- 
Jp^P^i^da' = (-1)1-1 ||^ + [^[ji /p^<|rfa (B.24) 



(B.25) 
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for ?7i < 

Finally, assembling the previous results, we have for the closed case that : 

2{l' + my. 2{l + m)\ 



•^ T 



{21' + 1){1' - my. {21 + l){l - my. 



^"^^^n-^^rifc.^.iMM'.m) 



4 

9=0 



(P)c,*(P) 



Srnm'^ / d^;V(2L+l) , I"' "^\'^ V {k, w)e{lj' , m) 
2 J ^ \b{k,w)Y 



(B.26) 



SP) 



where we have used the functional form (|147p for the function c^j . Recall that 
k^ — KL{L + 1) in the closed case. The function e{l, I', m) reflects the normalization 
of the spherical harmonics and allows for the particular case m = 0. It is defined as 
follows: 



.r//'^->m- i'2i + m2i' + m-my{r-my 



e(M',^^0)^4,/&lK^ (B.28) 



eU,i,m<Uj-y ^g^2 [{l + my.{l' + my) ^^'^^^ 

In the open case we have that (where k^ = \K\{X^ + j}): 

1,,^, /", /",,,, , / x^c 2{r + my 2{l + my 

< ail m' aim > = 7^1-f^l / "w; / dAA tanh (tt A) d; 



27r' 'J J ' ' """ {2U + 1){1' - my {21 + 1){1 - my 

\m\ 

■7r^Y[4Fi\P{k,w,t)e{l,l',m) 

g=0 

= 6mm'2TT\K\ / dw / dAAtanh(7rA) .I"' "'\' P{k,w)e{l,l' ,m) 
J J \d{k,w)\^ 

(B.30) 
Finally in the flat case we have that: 

r r Q(P)Q*iP) 

< ai'rn'ai^ > = 5^,n'2TT / dw / kdk-ff- — ^V{k,w)e{l,l',m) 
J J \d{k,wy^ 

(B.31) 
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Appendix C. Gauge transformations 

We assume that we are working with theories derived from generaUy covariant actions. 
Consequently, we may use the four dimensional diffeomorphism invariance to impose 
restrictions on the pcrturbational quantities in order to simplify calculations. Consider 
an infinitesimal vector field S,^, which will be used to generate diffeomorphisms. In 
this sense infinitesimal means that all components of ^^ are to be regarded as first 
order in smallness. We may define its components as follows: 

^^dxf" = a'^'^i^tSdt + £,,£dz + ^pWaSdx" + U^aSdx") (C.l) 

k.m 

Under an infinitesimal diffeomorphism the metric changes at fixed coordinate 
location changes as follows: 

gt^u -^ g,,u + ^D^^^ (C.2) 

where D^ is the derivate operator compatible with the background spacetime 
metric. 

We now decompose the tensor D(^^S,u) into polar and axial parts which are, 
respectively 



+ 2 (^i + dz^t)£dtdz + -{d,^p + ^,)Vb£dzdx' 



bouzax 

^pVaybSdx'^dx'' - n^aab£dx''dx'' 

(dz^z - n^t)£dzdz) (C.3) 



{D^^^)^'^Ux''dx'' = a^ J^'yl^-^A'^bSdzdx'' + ^(^'^S/bSdtdx'' 



'2 ' 2- 

k.m 



+ U^iayb)£dx''dx'') (C.4) 

Recall that the perturbed metric has the following polar and axial components: 



Sglf^Ux^^dx" = a^ ^{V{z, t)£dt^ + E{z, t)£dzdt + F(z, t)£d. 



z" 



k,m 



+ B{z, t)Vd£dzdx'' + C(z, t)Vd£dtdx'' 

+ U{z, t)-iab£dx''dx^ + X{z, t)\Iayb£dx''dx^) (C.5) 

and 

SgliUx'^dx'' = a^ ^{R{z,t)Va£dtdx'' + S{z,t)V a£dzdx°' 



k.m 



2Q{z,t)V(aVb)£dx''dx'') (C.6) 
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It is clear from the form of (jC.31) and (|C.4I) that none of the individual components 
in (|C.5I) and (IC.6P are unalterable under a gauge transformation. However, we may 
construct combinations of perturbations that are gauge invariant. Combinations in 
the polar case are: 

Hr>=»i.-| + f + ^ (0.7) 

Ef^ =nE-nB' + iwnC + iwU (C.9) 

4'^-B'-':^-E-^:^ (c.io) 

Combinations in the axial case are: 

E[^^^iwQ-S (C.ll) 

4^^^ == iwi? - S" (C.12) 

Appendix D. Comparison to FRW scalar-tensor- vector decomposition in 
flat limit 

The shearless class of spacetimes described by © allow only spatially flat FRW 
spacetimes as a limiting {K — )• 0) case. 

In a general gauge we may write a perturbed spatially flat FRW spacetime as 
follows: 



ds^ = a2(-(l + 2*)df2 + (v,a + a,)dtdx' + (% + 5g.,j)dx'dx^lD.l) 

6gij = - 2$% + ViVj/3 + V(,;3j) + /ly (D.2) 

where Via* = Vi/3' = ft.'j = V^/i* = 0. Indices i and j are three-dimensional, 
indices are raised and lowered with the background co- moving spatial metric Sij, 
and Vi = di is the metric compatible derivative operator associated with Sij. The 
perturbations *,$,a,/3 are termed scalar perturbations, the divergenceless fields /3* 
and a' are termed vector modes each with two polarizations, whilst the field hij carries 
two polarizations. 

We now write a,; and /3i in terms of the 2+1 polar and axial decomposition: 

a. = ^a.£' (D.3) 

(D.4) 
(D.5) 
(D.6) 





k,m 


aa -- 


= Y^iCVaS + C,Va£) 




k,m 


Pz'- 


k,m 


/3a = 


k,m 
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where recall that in the flat case £" is a Bessel function J,n(kr). 

The vanishing divergence of a^ and /3j then yields the following conditions: 

= V,a, - k\ (D.7) 

= V^/3, - k^Tj (D.8) 

Similarly we can decompose the field hij as follows: 



hijdx'^dx^ = ^ hzzS dzdz + h^adzdx"' + habdx'^dx^ (D.9) 

hza - ^(/i^^Vaf + /i(^)v„£:) (D.IO) 

k,7n 

Kb = ^(/l(^)7afc + /^^''^VaVbf + /l(^)V(,Vfc)£:) (D.ll) 
k.m 

The condition /i*j = then implies that: 

= hzz + 2h^^^ -k^h^^^ (D.12) 
whilst the transverse condition 9i/i' — implies that: 

= dzh\ - k^h^^'^ (D.13) 

= (a,/i(^)-fc2/i(^))Va£: (D.14) 

= (a,/./^) - y/^(^')V.£- (D.15) 

We would like to see how the polar and axial perturbations are related to tensor, 

vector, and scalar perturbations in the limit of spatial fiatness. Comparison between 
(|D.1|) and (jC.Sp and (|C.6p yields, in terms of Fourier modes: 

y = -2* (D.16) 

E = iwa + az (D.17) 

F = - 2$ - u;2/3 + iwl3z + h^^ (D.18) 

B =m;/3 + -(iw77 + /?^) + ft.(^) (D.19) 

C=a + C (D.20) 

f/ = h^T) _ 2$ (2 21) 

X = /3 + 77 + /i(-°) (D.22) 

and 



(D.23) 
(D.24) 
(D.25) 



R 


= c 


S 


= —iwrj + Iv- ' 


2Q 


^^r^ + U^) 
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We see then that in the flat hniit, axial components consist entirely of tensor and 
vector perturbations whilst polar components are themselves sums of scalar, vector, 
and tensor perturbations. 

Consider the gauge invariant combinations of the previous section. We see that 
the following two are composed entirely of tensor perturbations in the flat limit: 

eP ^ ^wh^^^ -hll + ^+^^^ (D.26) 

1 2 2 2 ^ ' 

e[^) ^ ^yjl^ _ hiV) (D.27) 
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Appendix E. Table of selected notation 



Symbol 



Definition 



"fab 

flij 

\K\ 
t 
z 

x" 
ait) 
bit) 
D, 

Va 
V2 

SG 






5G' 

ck 
<^m 

W 

<Piz,x\t) 

e 

r 

a 

pa 

A 

Tit) 

Q{P)ix\P\t) 

e^'^^x^P^t) 

'-'na 

'-'na 
Q^Pjk,w,t) 

e^j^lik,w,t) 

^NA 
fiab..M 

Qab..M 

o(p) 

(A+)° 

\yab..M 



Z' 



ab..M 



cd..M 
cd..M 



Co-moving 2-surface background metric 

Co-moving 2-surface background volume form 

Background volume element on the co-moving 2-surface 

Background Riemann curvature tensor on the co-moving 2-surface 

Background 3-dimensional spatial metric 

Co-moving background 3-dimensional spatial metric 

Curvature scalar of the co-moving 2-surface background metric 

Conformal time 

Co-moving coordinate in the z-direction 

Comoving coordinate on the co-moving 2-surface 

The background expansion rate of the spatial surface orthogonal to idz)'^ 

Background expansion rate along the z-direction 

Derivative operator compatible with background spacetime metric 

Derivate operator compatible with co-moving 2-surface background metric 

Laplace-Beltrami on the background co-moving 2-surface 

Component of polar perturbation to the spacetime metric 

Component of axial perturbation to the spacetime metric 

Component of the polar perturbation to the Einstein tensor 

Component of the polar perturbation to the Einstein tensor 

Eigenfunction of V^ corresponding to eigenvalue —kk* for generic curvature 

Wavenumber in the z-direction in the shearless case 

Defined as f-a^b 

Scalar field supporting the shear-free anisotropic curvature 

Infinitesimal difFeomorphism generating vector field 

Distribution function for species of massless particles ^ 

Distribution function for species of massive particles 

Four momentum vector of particle (massless or massive) 

The norm-squared of the particle spatial momentum 

The projection of the co-moving spatial momentum along idz)^ 

Components of momenta in the 2-surface with normalization 'yabP°"P^ — 1 

AfRne parameter of null-trajectories 

The background photon temperature 

The dimensionless polar temperature perturbation 

The dimensionless axial temperature perturbation 

Basis function for polar mode-moment expansion 

Basis function for axial mode-moment expansion 

Polar temperature perturbation moment 

Axial temperature perturbation moment 

Normalization function for polar temperature moment 

Normalization function for axial temperature moment 

Polar Chebyshev tensor (M indices) 

Axial Chebyshev tensor (M indices) 

Function appearing in the Boltzmann equation 

Matrix relating directional derivatives of Bna to sums of such functions 

Polar projection operator 

Axial projection operator 

2- vector defined via equation ([57]) 

2- vector defined via equation (|88l) 

Tensor defined via equation (|78l) 

Tensor defined via equation (I120p 
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Symbol 


Definition 


/? 


Cosine of angle between z direction and 3-d wave-vector in the spatially flat limit 


Ja{x) 


Spherical Bessel function 


T 


Particle proper time 




Spherical harmonic 


i(+)a 


The equivalent of equation ()87l) but instead defined with respect to ^cd 


i(-)a 


The equivalent of equation (|88p but instead defined with respect to jcd 


rpfyn) 


Matrix element for a representation rj of a group 


V{k,w,t) 


Dark matter power spectrum 


■pii) 

1 mn 


Function defined via eauation (IA.7|) 


prn 


Associated Legendre function 


Va 


Legendre polynomial of order A 


'D+ 


Differential raising operator fsee for instance equation (lA. 551)1 


2?_ 


Differential lowering operator fsee for instance eauation (|A.56I)') 




Eigenfunction of the co-moving 2-surface Laplace-Beltrami in the open case 


Function defined via eauation (|A.4ip 


^n—m\p^) 


Bessel function of the first kind 


jk 


Function defined via eauation (|A.70I) 


aim 


Polar temperature perturbation in spherical harmonic space 


S(P) 


Gauge-invariant polar perturbation 


^(A) 


Gauge-invariant axial perturbation 



